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SYNOPSIS 


The most labour intensive, uncontrollable area in the manufacturing of precision 
parts involves machining operations, which frequently demand as much as 15 % of the 
total manufacturing cost. The result of high quality finish on the parts is improved 
performance, running-in-period’ of minimum duration, and considerable increase in the 
length of life of the component. Abrasive finishing process is usually employed as the last 
operation on the high value added components. It is imperative to achieve these processes 
economically Ynth strategic consideration of the key influencing parameters. 

Grnnding, lapping, honing, and super finishing are traditional methods of finishing. 
But the applications of grinding, lapping, honing, super finishing, and other abrasive 
finishing processes are limited to the production of the workpieces of the basic forms such 
as flat, cylindrical etc. These conventional finishing processes are being pushed to their 
limits of performance and productivity in general, and in case of hard materials, and 
components of complicated shapes in particular. Hence, a need is being felt to develop a 
finishing process with wider bounds of application areas, better quality performance, and 
higher productivity. Abrasive flow machining is one of such processes. 

Abrasive flow machining (AFM) is a non-traditional finishing process that is used 
for deburring, radiusing, surface finish improvement, polishing of irregular shaped passages, 
and removal of thermal machining recast layer. With today’s focus on total automation in 
the flexible manufacturing system, AFM process offers both automation and flexibility in 
final machining operations. A key feature differentiating AFM from most of the other 
finishing processes is the ability to control and select the intensity and location of abrading 
through fixture design and abrasive media selection. Applications of this process include 
finishing of components in aerospace, automotive, dies, chemical processing and medical 
industries. Materials from soft aluminum to tough nickel alloys, ceramics and carbides have 
been successfixlly micro machined by this process. AFM process provides a high level of 
surface finish and close geometric tolerances with economically acceptable rate of surface 
generation for a wide range of industrial components. The ability of media in AFM process 
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to finish diflScult to reach areas, to follow complex contours and to simultaneously work on 
multiple edges and surfaces, makes it more versatile than other finishing processes. 

AFM removes small quantity of material by flowing a semisolid abrasive laden 
compound called 'media* (abrasive particles uniformly suspended in viscous chemical 
compound) through or across the surfaces of the workpiece to be finished. Two vertically 
opposed cylinders extrude abrasive media back and forth through passages formed by the 
workpiece and tooling. The machining action compares to grinding or lapping operation as 
the media gently and uniformly abrades the surface or edges. The media acts as a 'self 
deformable stone' having protruding abrasive particles acting as cutting tools. 

AFM process is very much in its early stages of development. Sustained research is 
required to transform the process into a mature technology, for its successful induction into 
the industry. The thesis makes an attempt in this direction. 

Both theoretical and experimental studies of abrasive flow machining are greatly 
hampered by the inherent random nature and multiplicity of variables. Not much 
information is available in the literature which deals with theoretical analysis of material 
removal mechanism in AFM. The relationship between process parameters and 
performance characteristics of the process are not completely known. There is a need to 
carry out systematic theoretical analysis of viscous media flow to determine the stresses and 
machining forces acting on the workpiece. This will help in modeling and proper 
understanding of the mechanics of the micro cutting process during AFM. The objective of 
the present work is to investigate the mechanism of AFM process and to study the effect of 

various process parameters on the process performance by experimentation and theoretical 
modeling. 

The specific objectives to achieve the above mentioned goals are as follows : 

• To conduct experiments to study the effects of various process parameters on material 
removal and surface roughness. 

• To analyze the flow of media in AFM process using finite element method by 
considering the media as viscoelastic material. 

• To develop a model to determine the material removal and surface finish in AFM 
process and to validate the developed model with the data obtained by experimentation. 
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• To develop a model for the topography of media and estimate the active grain density 
by microscopic technique; and to simulate the random surface generated in AFM 
process. 

• To analyze the specific energy and thermal aspects of AFM process. 

• Application of neural network in AFM process for modeling and optimization. 

Experiments have been conducted on the AFM machine designed and fabricated at 
Indian Institute of Technology, Kanpur. Experiments have been planned using central 
composite rotatable design to get the useful inferences by performing minimum number of 
experiments. A fixture was designed to direct the media to flow in the desired areas. The 
process parameters considered are number of cycles, extrusion pressure, percentage 
concentration of abrasives and reduction ratio. Internal surfaces of mild steel workpieces 
with circular cross section have been finished using abrasive media (SiC abrasives mixed 
with polyborosiloxane carrier). Material removal increases with extrusion pressure, 
concentration of abrasives, reduction ratio and number of cycles. Surface roughness value 
decreases with increase in extrusion pressure, reduction ratio and concentration of 
abrasives for a specified number of cycles. 

A model has been proposed for analysing media flow through the media cylinder 
and workpiece. The method involves the use of finite element technique to solve the 
governing equations in terms of primary variables (velocity and pressure). The theoretical 
analysis of media flow predicts the stresses developed during the flow. These computed 
stresses can be used to evaluate the machining forces. The model predicts the normal 
stresses at the workpiece surface with reasonable accuracy. 

AFM is considered to be a scratching action performed by the abrasive grains in the 
flowing media. A model has been proposed to determine the material removal and 
generated surface roughness during AFM using abrasion theory. The model is based on 
accumulated plastic flow of material by repeated indentation by moving abrasive particles. 
The model includes the operating conditions of the process, e.g. abrasive size, extrusion 
pressure, workpiece material hardness and size of the workpiece to be machined. A good 
agreement is found between the experimental results and theoretical predictions. Material 
removal increases with increase in extrusion pressure and percentage concentration of 
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abrasives in the media. Surface roughness value decreases with increase in extrusion 
velocity / extrusion pressure, percentage and concentration of abrasives for a specified 
number of cycles. 

Understanding of the process of generating finished surface is important for the 
control of the AFM process. Systematic modeling and analysis of AFM are greatly 
hampered by the inherent random nature and multiplicity of variables. Simulation technique 
can help to understand the mechanics of AFM process by allowing the effects of various 
process parameters. The concept of a newly developed, versatile simulation model which 
generates and statistically evaluates the interaction between spherical abrasives and 
workpiece surface is described. The simulation enables prediction of the roughness of the 
machined surfaces and material removal with reference to percentage concentration of 
abrasives, extrusion pressure of media, reduction ratio and number of cycles. The number 
of active grains and state of cutting edges are closely related to the AFM performance. The 
present stochastic simulation method is used to evaluate the active grain density. A 
microscopic technique is also used to estimate the number of active grain density at various 
mesh sizes and concentrations. The results obtained from stochastic simulation method 
agree well with the data obtained by microscopic technique. 

A model has been proposed for the determination of specific energy and 
tangential forces in AFM process. The model includes the operating conditions of AFM 
e.g., grain size, applied pressure, hardness of workpiece material, number of cycles and 
active grain density. A simple model of heat transfer in AFM has also been developed 
that considers heat flow to the workpiece and abrasive media. The dependence of the 
workpiece temperature on the AFM parameters has been explored. The model also predicts 
the fraction of heat entering into the workpiece and media. 

AFM process is too complicated to warrant appropriate analytical models and most 
of the time, analytical models are based on many assumptions, which contradict in reality. 
More importantly, it is sometimes difficult to adjust the parameters of the above mentioned 
models according to actual situation of the machining process. Therefore, neural networks 
which can map the input/output relationships and possess massive parallel computing 
capability, have attracted attention of researchers working in the area of AFM process. 
The thesis presents architecture of back-propagation neural network model for 
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prediction of material removal rate and surface roughness in AFM. In order to decide the 
structure of neural network, the rate of error convergence is checked by changing the 
number of hidden layers and also by adjusting the learning rate- and momentum rate. The 
three layers back-propagation with four inputs, two outputs and nine hidden nodes is 
adopted for neural network. It is observed that learning could remarkably be enhanced by 
training the network with noise injected inputs. The results confirm the feasibility of this 
approach, and show a good agreement with experimental results for a wide range of 
machining conditions. 

The neural- network has also been used for optimum selection of input 
parameters of AFM process. First, a generalized back-propagation neural network is used 
to establish the process model. A second network, which parallelizes the Augmented 
Lagrange Multiplier algorithm, determines the corresponding optimal machining parameters 
by minimizing a performance index subject to appropriate operating constraints. An 
important consideration is that process optimization can be performed in the absence of 
process model and purely by experimental observations. 

To validate the neural network results, optimization of the AFM process has been 
also performed using genetic algorithm (GA). Genetic algorithms are computerized search 
and optimization algorithms, based on the evolutionary process of biological organisms in 
nature. The optimization results of neural network coincide well with the results obtained 
by genetic algorithms and justify the neural network approach. 

The thesis also outlines the possible avenues for future research in abrasive flow 
machining process. 

Organization of the thesis 

First chapter of the thesis comprises introduction of traditional finishing methods 
and their limitations, introduction and importance of abrasive flow machining process, 
review of literature, and scope and objective of the present work. Description of the 
experimental set-up, details of experimentation and the associated results constitute chapter 
2 of the thesis. Chapter 3 presents a model for analysing the flow of AFM media using 
finite element technique. Modeling for material removal and surface roughness using 
abrasion theory is also described. A model for the topography of abrasive media and 
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estimation of active grain density by microscopic technique is described in chapter 4. A 
modified model to simulate the random surface generated and material removal in AFM is 
also discussed. Chapter 5 is devoted to forces, specific energy and thermal aspects of AFM. 
Chapter 6 deals with the modeling of AFM by neural networks. Optimization of AFM 
process by neural network and genetic algorithm is also presented. The conclusions drawn 
from the present work are enumerated in chapter 7. Possible avenues for future research in 
AFM are also outlined. The references and appendices are given at the end of the thesis. 
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Chapter 1 

Introduction and Literature Review 


1.1 Abrasive Fine Finishing 

The most labour intensive, uncontrollable area in the manufacture of precision 
parts involves machining operations, which frequently demand as much as 15 % of 
the total manufacturing cost The cost of surface finish increases sharply for a roughness 
value of less than one micron. Surface finish has a vital influence on important functional 
properties such as wear resistance and power loss due to friction on most of the 
engineering components. When constant fiiction is required, the surfaces must be uniform 
and have a predictable wear rate. Surface finish can also affect part-to-part fit and 
tolerances. The result of high quality finish on the parts is improved performance, 
‘running-in-period’ of minimum duration and considerable increase in the length of life of 
the component. These final machining operations are uncontrollable and if they are not 
carried out properly, performance of the machined component may be far from optimum. 
The dimensional and alignment accuracy and the quality of the surface finish are taken 
care of by finishing processes. 

The need for high accuracy and high efficiency machining of difficult-to-machine 
materials is making the application of abrasive finishing technologies increasingly 
important Abrasive finishing process is usually employed as the last operation on the 
highly value added components. It is imperative to achieve these processes economically. 
This is generally possible with strategic consideration of the key influencing parameters. 

Traditionally, single or multiple point tools with geometrically defined edges, such 
as turning and milling tools, have been considered for roughing and semi-finishing. For 
fine finishing, fine abrasives (be they loose or fixed) are invariably used in the components 
of the highest quality in terms of form and finish accuracy, and surface integrity. The basic 
idea of abrasive fine finishing processes is to use a large number of multipoint or random 
cutting edges with indefinite orientation and geometry for effective removal of material in 
the form of chips having sizes smaller than those obtained during finishing methods that 
use cutting tools with defined edges. Because of extremely thin chips produced in abrasive 
machining, it yields better surface finish, closer tolerances, more localized control. 
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generation of more intricate surface features, and machining of harder and difFicult-to- 
machine materials [1,2] Unlike other machining processes, the failure of one cutting edge 
does not affect the process in abrasive machining. The abrasive tools are consequently 
quite stable. Compared with the defined cutting edge tool, the abrasive tool is less 
sensitive to interrupted cuts and variations in work hardness as well as the areas of the 
workpieces to be machined 

Grinding, lapping, honing, and super finishing are commonly used methods of 
finishing. Lately, magnetic abrasive finishing (MAF) has also joined this class of abrasive 
finishing processes. Even though these processes are in use for various applications, each 
process has its limitations in producing the desired surface finish on the components. 
Some of them are discussed in the forth-coming sections. 

1.2 Abrasive finishing processes 

Grinding [3] is a finishing process, which uses abrasive wheel or tool made of 
large abrasive grains retained by bonding material. Although grinding is more efficient for 
removing material than other finishing methods, it is still difficult to achieve desired fine 
finish by grinding alone. Though the finish can be improved with the application of 
grinding wheels with fine grits, they get heavily loaded with debris during the grinding 
process. Hence, very fi-equent dressing is needed to remove the loading, which causes 
excessive wheel loss and interruption in production. Finishing of intricate parts by 
grinding requires expensive profile grinding wheels. Grinding of certain materials such as 
aluminum and brass is not possible as they cause quick loading of the grinding wheel, 
making grinding quite difficult. Grinding develops localized heat with potentially harmful 
effects on the integrity of the surface. 

Lapping [4] is a low pressure abrading process, which is used for obtaining 
surface finish, and dimensional accuracy superior to that obtained by grinding and honing. 
Relatively low force is exerted on the workpiece that is usually held in a well supported 
manner, consequently lapping is often successfully applicable to brittle materials and 
fi-agile parts. The rate of stock removal is low due to low cutting speed and shallow 
penetration of fine abrasive grains into the workpiece surface. The material is removed by 
applying loose abrasive slurry between work material and close fitting surface called lap 
(Fig. 1.1) The abrasive grains follow a free continuously changing path over the 
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workpiece and produce a stress free surface by cutting action [5], The shapes of the 
surfaces frequently worked by lapping are mostly limited to basic forms, such as flat, 
cylindrical (including slender tapers), and spherical. In lapping, to polish a three 
dimensional surface such as finishing of dies, the operation is usually manual. Also the 
accurate motions of lap and work are required, necessitating use of high precision set-up 
and linear actuators, increasing the cost 



Grain 2 


Fig. 1.1 Schematic diagram of a lapping process 

Honing [6,7] is a low velocity abrasive machining process characterized by low 
cutting pressure and large area of contact The principle of honing, with the self- 
piloting tool automatically finding the neutral axis of the hole and developing the new 
surface in relation to that axis, are well-adapted for producing straightness and roundness 
by removing the least amount of material often within very short cycle time. Bonded 
abrasive stones are forced against the work surface, while a combination of rotary and 
traversing motions successively and repeatedly bring the entire work surface into contact 
with the stones (Fig. 1 .2). It is employed to improve the accuracy of cylindrical, usually, 
internal surfaces. It also reduces or corrects geometric errors resulting fi-om previous 
operations, such as machining, heat treating etc., with least possible amount of material 
removal. It is used in a wide variety of applications including engine cylinder walls, 
compressor bodies, valve bodies, bearings and hydraulic cylinders. The surface 
topography generated by honing can have a profound effect on tribological performance 

m. 
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Superfinishing [4] is an operation using bonded abrasive stones in a particular 
way to produce an extremely high quality surface finish. A very thin layer of material is 
removed in superfinishing. In this process, a very fine abrasive stick is fixed in a suitable 
holder and pressed against the surface of the workpiece with a light spring pressure. The 
stick is given a feeding and oscillating motion, and the workpiece is given rotating or 
reciprocating motion according to the requirement of the shape being machined In 
superfinishing operation, there occurs formation of sharp edge which is subjected to 
manual working. 

In magnetic abrasive finishing process [ 8 ], granular magnetic abrasives 
composed of ferromagnetic materials (as iron particles) and abrasive grains (as AI 2 O 3 , SiC 
or diamond) are used as cutting tools and the necessary finishing pressure is applied by 
electro-mechanically generated magnetic field (Fig 1 3). The magnetic abrasives are 
joined to each other magnetically between magnetic poles along the lines of magnetic 
force, forming flexible magnetic abrasive brushes. When a workpiece with rotatory, 
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vibratory, and axial movement is inserted in such a magnetic field, surface and edge 
finishing are performed by these magnetic abrasive brushes. MAF can be used for 
finishing of tubes, external finishing of rods and finishing of flat surfaces made of 
ferromagnetic and non-ferromagnetic materials. 



Fig. 1.3 Schematic view of magnetic abrasive finishing (MAF) 1^1 

But the applications of the grinding, lapping, honing, super finishing, magnetic 
abrasive machining and other abrasive finishing processes are limited to the production of 
the workpieces of the basic forms such as flat, cylindrical etc. These finishing processes 
are being pushed to their limits of performance and productivity in general, and in case of 
hard materials and components of complicated shapes in particular. Hence, a need is 
being felt to develop a finishing process with wider bounds of application areas, better 
quality performance and higher productivity. Abrasive flow machining is one of such 
processes. The characteristics of four important finishing processes are compared in 
Table 1.1. 

The foregoing discussion clearly indicates that abrasive machining has a vital role 
to play in the finishing of industrial components. The work reported in this thesis is on a 
new process, called abrasive flow machining (AFM). The following section provides an 
introduction to this process, and a review of the relevant literature available. 
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Table 1.1 Comparison of abrasive finishing processes. 


s. 

N 

Process 

features 

Lapping 

Honing 

MAP 

AFM 

■ 




0 04 [8] 

0 05 [17] 

1 

Dimensional 

tolerance, 

(//m) 

±0.5 [1] 

±0.5-1.25 [1] 

±0.5 [8] 

±5 0(17] 

3. 

Material 

removal, (mm) 

<0 0025[1] 

0.061-0.183 [1] 

0.002 - 0 007 [8] 

0.008-0.010 [9] 

4. 

Pressure, 

(MPa) 

0.01-0.2 [1] 

1.0 -3.0 [1] 

0-0.077 [8] 

0.69 - 22.0 [9] 

5. 

Abrasive 

product 

type 

Abrasive grain 

entrained in a 

liquid vehicle 

Bonded 

abrasives 

1 

Magnetic 

abrasives 

Composed of 

ferromagnetic 

particles and 

Conventional 

abrasive grits 

Semisolid 

abrasive media 

composed 

of viscoelatic 

carrier and 

abrasive grits 

1 

Work surface 

Configuration 

Flat, 

cylindrical, 

and spherical 

surfaces 

Cylindrical i 

Surfaces i 

Flat and 

cylindrical 

surfaces 

Inaccessible areas 

and complex 

Internal 

passages 


1.3 Abrasive Flow Machining (AFM) Process 

With today’s focus on total automation in the flexible manufacturing system, the 
abrasive flow machining process offers both automation and flexibility in final machining 
operations. AFM can be applied to an impressive range of finishing operations, providing 
uniform repeatable and predictable results A key feature differentiating AFM fi-om most 
finishing processes is the ability to control and select the intensity and location of abrading 
through fixture design, media selection and process parameters. Abrasive flow machining 
is a non-traditional finishing process that is used to attain one or more of the following 
objectives: 

• Uniform radius generation. 
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• Edge or intersection deburring 

• Surface finish improvement 

• Polishing of irregular shaped passages 

• Removal of thermal machining recast layers 

• Post processing of conventionally machined surfaces to eliminate residual stresses. 


Applications of this process include finishing of components in aerospace, 
automotive, dies, chemical processing and medical industries. It is an appropriate process 
for pharmaceutical and semiconductor processing industries. Materials fi'om soft 
aluminum to tough nickel alloys, ceramics and carbides have been successfully micro 
machined with the process [9]. AFM process provides a high level of surface finish and 
close geometric tolerances with economically acceptable rate of surface generation for a 
wide range of industrial components. The ability of media in AFM process to finish 
difficult to reach areas, to follow complex contours and to simultaneously work on 
multiple edges and surfaces, makes it more versatile than other finishing processes. 
Abrasive flow machining is a cool process, as opposed to grinding which develops 
localized heat with potentially harmful effects on the integrity of the surface. 


Fixture 



Upper cylinder 

Abrasive media 

Tcx5ling 

Workpiece 


Lower cylinder 
Piston 


Fig. 1.4 Schematic diagram of abrasive flow machining process 
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AFM removes small quantity of material by flowing a semisolid abrasive laden 
compound called 'media* (abrasive particles uniformly suspended in viscous chemical 
compound) through or across the surfaces of the workpiece to be finished [10]. Two 
vertically opposed cylinders extrude media back and forth through passages formed by the 
workpiece and tooling (Fig. 1.4). The machining action compares to a grinding or lapping 
operation as the media gently and uniformly abrades the surfaces or edges. The media acts 
as a 'self deformable stone' having protruding abrasive particles acting as cutting tools. 
The media is composed of semisolid carrier and abrasive grains. AFM systems designed 
for production applications often include part cleaning and unload/reload stations, media 
maintenance devices and cooling units. 

In order to achieve the desired machining, the following process-influencing 
parameters should be considered. 

1 . Setting parameters: Extrusion pressure, velocity of flow of media, flow volume, and 
number of cycles. 

2. Tooling : Shape and size, stiffness, stability and resulting media flow passage. 

3. Media : Type and size of grains, concentration of abrasives, and type, rheology and 
viscosity of carrier. 

4. Workpiece : Hardness, surface roughness, shape and stiffness. 

The most commonly used abrasive grains are aluminum oxide (AL2O3), silicon 
carbide (SiC), boron carbide (B4C) and diamond. Sometimes their mixture can also be 
used. The carrier is made by blending rubber like polymer with a gel like lubricating 
diluent, which can give desired velocity e.g. polyborosiloxane. Tooling design directs the 
media to flow in the desired area. 

1.3.1 Unique characteristics of AFM process 

AFM offers a number of unique and important features. Some of them are reported 
below [9,10,11]. 

1 . Deburring, radiusing and polishing in one operation even in inaccessible areas. 

2. Produces true round radii even on complex edges. Since the media confirms to the 
passage geometry, complex shapes can be processed as easily as simple one. 

3. Reduces surface roughness by 75 to 90 percent on cast, machined or EDM'd surfaces. 
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4. Confidently produces close tolerances due to the selectivity and easy controllability of 
the process. 

5 Processes dozens of holes or parts simultaneously in about the same time it would take 
to do just one. 

6 Multiple passages can be processed simultaneously with uniform results. Cooling air 
holes on turbine disks and hundreds of holes in a combustion liner are deburred and 
radiused in one operation. 

7. APM maintains flexibility. The same machine can do a variety of jobs just by 
changing tooling, process settings and, if necessary, abrasive media. 

8. Defects are removed or blended in pure grinding mode. There is no peening over or 
smearing of metal into the surface. 

9. Jobs requiring hours of highly skilled hand polishing are processed in a few minutes 
automatically. 

10. Process variables are easily controlled. Small changes in media flow rates, pressures, 
or volumes normally have little effect on results. 

1.3.2 Applications of AFM 

AFM is used for finishing of various components in wide range of industries. 

Some of them are given below [9, 10, 1 1]. 

(a) Aerospace 

• Deburring of aircraft valve bodies and spools. 

• Polishing and radiusing of blades, disks, hubs, shafts and air cooling holes in the 
turbine disk. 

• Finishing accessory parts such as fuel spray nozzles, fuel control bodies and bearing 
components. 

• Removing the recast layers in laser'd and EDM'd cooling holes on turbine disk and 
blades. 

(b) Dies 

• AFM has been successfully applied to polish extrusion, drawing, wire-forming, 
compacting and cold-heading dies. 

• The process polishes surface in the direction of flow, producing a better quality and 
longer lasting dies with uniform surface and gently radiused edges. 
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(c) Automotive 

• For deburring and polishing fuel injection systems, splines, gears, pumps, valves and 
fittings. 

• AFM abrades, smoothes and polishes the surface of the two stroke cylinders and four 
stroke engine heads for improved air flow and better performance. 

(d) Chemical processing 

• To finish valves and fittings used for pharmaceutical and semi-conductor applications. 
AFM can easily reach inside the smaller fitting and valves, where access for hand 
polishing is very difficult. 

Studley [12] presented applications of AFM process, e g. finishing of hundreds of 
holes in a disk in one operation, uniform and controlled stock removal of cast blades with 
minimum dimensional change, resizing of cooling holes in aircraft engine combustion 
liners, finishing of turbine components which otherwise are very difficult to finish. The 
AFM process has been successfully applied for the components ranging from gears as 
small as 1.5 mm in diameter and orifice as small as 0,2 mm, to splined die having passages 
of 50 mm in diameter [9], Some of the applications of AFM process for finishing various 
components in different industries are shown in Fig. 1 ,5 to F'ig. 1.8. 



Fig. 1.5 Beburring and polishing gear teeth [lO] 
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Fig. 1.6 polishing of airfoil and gas path surfaces in turbine engine parts [15] 



Fig. 1.7 Polishing of both tool steel and tungsten carbide dies [15] 



Fig, 1.8 Valves and fittings used for pharmaceutical and semiconductor 
applications [15] 
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1.4 Literature Review 

The Abrasive Flow Machining (AFM) process has been developed and brought to 
production stage in U. S A. during the years of 1966 to 1968. The process was initially 
used to perform critical deburring of aircraft valve bodies and spools, and production of 
burr free internal edges In 1978, AFM was reported as an alternative for hand polishing of 
extrusion dies at the International aluminum extrusion technology seminar in Atlanta, 
Georgia. Since that time, use of AFM for polishing of dies has developed as a way of 
satisfying stringent manufacturing requirements. AFM process is in its early stage of 
development; and the available published literature mainly deals with a qualitative 
description of the process and its applications. The review as given below, is presented in 
two parts, viz related to experimental findings and theoretical investigations. 

1.4.1 Review of experimental studies 

Many researchers conducted experimental studies to investigate the effects of the 
process parameters on the process performance. The effects of process variables like 
extrusion pressure, number of cycles, viscosity, concentration and mesh size of the 
abrasive media are discussed 

Rhoades [9,10,11] studied the basic principle of AFM process and identified its 
process control parameters. The depth of cut primarily depends upon the size, relative 
hardness, sharpness of abrasive grains and extrusion pressure. The number of cuts made 
by abrasive particles located on the outer periphery of the slug (the media extruded 
through the workpiece, shown in Fig. 1 .4) during one cycle, is a function of length of the 
media as it passes through the restriction, and size and concentration of abrasives in the 
media. The type of flow pattern to occur has been shown to depend upon the machine 
setting, media formulation, workpiece material and tooling configuration. 

It is also reported [10,11] that higher viscosity media is to be used for uniformly 
abrading the walls of large passages and lower viscosity media is to be used for radiusing 
edges and for processing small passages. The effective life of the media depends on a 
number of factors including the initial batch quantity, the abrasive size and type, the flow 
speed and part properties and its configuration. 
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Przykienk [13] described the effects of workpiece geometry, surface roughness and 
different media on process performance. He concluded that with small bore diameter more 
grains come in contact with the wall and material removal increases. The edge radius 
decreases if diameter and length of a bore increase. The share of abrasives in media, the 
grain size and viscosity of base media are the most important parameters which have 
influence on stock removal and media velocity. In comparison to this the influence of 
material, temperature, tjqpe of abrasive and bore geometry is small within the given limits 
It was also found that the media should not be allowed to achieve temperature more than 
100°C to keep the viscosity of the media within certmn limits. It is possible to use 
aluminum or polyurethane as material for the construction of fixtures. 

For any working pressure the amount of abrasion that occurs is directly related to 
the slug length of flow. It has been shown that if the two passages of different areas are 
given the same volume flow, then the smaller passage abrades more than the large passage 
due to greater slug length of flow [14]. 

Perry [15] presented some principles and practices regarding AFM. He reported 
that abrasion is high where media velocity is high. If two holes have different diameters 
but equal total flow area, then they will get abraded equally. High viscosity media 
provides high material removal rate. Low viscosity is used for radiusing purpose. Likewise 
coarse abrasives are used for faster cutting and fine abrasives for flowing through small 
holes, or for highly polished finishes. 

An increase in pressure and media viscosity increases material removal rate while 
surface finish value decreases [9,16-18]. Media viscosity combined with extrusion 
pressure relates directly to the media flow rate which affects the amount of abrasion, the 
uniformity of stock removal and edge radius size. Slow media flow rates are good for 
uniformly removing materials; high flow rates produce large radii [10]. It is found that for 
constant flow of volume, changes in media flow rate do not significantly affect material 
removal. With increasing number of cycles, material removal initially increases due to 
higher initial coarseness of the workpiece surface; but later on it decreases [17,19]. A 
study of the workpiece surfaces generated by AFM, with the help of scanning electron 
microscope photographs, indicates that there is significant improvement in the surfece 
finish only within initial few cycles [16]. Among all studied process parameters, the 
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dominating one is percentage concentration of abrasives in media followed by abrasive 
mesh size, number of cycles and media flow speed [19], 

The type of machining process used to prepare the specimens prior to AFM are 
found to significantly affect the improvement in surface finish [18]. The amount of 
material removed jfrom the WEDM’d and milled surfaces are significantly different from 
that of turning and grinding [18], because these machining processes produce different 
micro surface contours. WEDM’d surfaces are better suited to AFM than turned and 
milled surfaces. The actual reason for this suitability needs further investigation. However, 
suspicion goes to the weak junctions between the recast layer and the parent material 
formed during WEDM. 

Davies [20] reported the relationships between rheological behavior of 
polyborosiloxane-grit mixtures and their machining characteristics. It is observed from the 
experimental results that the relationship between number of extrusion cycles and changes 
in both the temperature and pressure drop depends upon the media type and 
polyborosiloxane-grit ratio. It is found that increase in the temperature results, decrease in 
viscosity of the media and increase in the volumetric flow rate. It has been observed from 
experimentation that viscosity of media increases with percentage concentration and 
decreases with the mesh size and temperature [21]. 

The physical properties of the media such as thermal conductivity, specific heat 
capacity and surface heat transfer coefficient of media have been experimentally 
determined by Fletcher [22]. These experimental values are also validated with the 
mathematical model of each property, which are then used in the calculation of thermal 
behavior of nuxture as a function of material composition and temperature. 

1.4.2 Review of theoretical investigations 

AFM process is complex due to little understood behavior of the media, and the 
complicated and random nature of mechanical action of material removal. William and 
Rajurkar [17] used a stochastic modeling and analysis technique called data dependent 
system (DDS) to study the AFM generated surfaces. The DDS modeling procedure 
consists of fitting ARMA models by standard non-linear least square routines until the 
residual sum of squares of random disturbances can not be significantly reduced as 
determined by F-test criterion. Wave length decomposition of AFM surface profiles has 
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indicated the presence of two wavelengths, a large wavelength linked to the main path of 
the abrasive grains and a small wavelength linked to the cutting edges. Methods have been 
proposed to estimate the number of dynamic active grains involved in the cutting and the 
amount of abrasive grain wear per stroke. The number of active grains has been obtained 
by multiplying the frequency associated with the primary root from the DDS models by 
the time per stroke and dividing by the cross sectional area of the extrusion passage. The 
range of active grains found is very close to 280 to 600 reported for the grinding 
operations. 

DDS modeling of the surface roughness profile data by Loveless et al. [18] 
indicated a general trend of AFM reducing the order of the ARMA models representing 
the surface. Green's function plots showed that AFM removed the signs of the original 
machining process of the workpiece surface. 

Rajeshwar et al. [23] presented a mathematical model to determine the 
characteristics of the media flow during machining. The simulation model has been 
developed using the constitutive equations of Maxwell model describing non-Newtonian 
flow characteristics of the AFM media through cylindrical workpieces. Finite difference 
method was used to obtain the solution. It is observed that there exists linear relationship 
between the shear stress on the surface and thickness of the material removed determined 
experimentally. 

Control systems may be added to monitor and control additional process 
parameters, such as media temperature, viscosity and flow speed of media, and wear of 
abrasives. The estimate of number of dynamic active grains may be useful as abrasive tool 
life criterion. It might also serve as an indicator that media has absorbed too much 
workpiece material and must be replaced [16]. These methods may aid in the off-line 
optimization of the process. They may also be useful as a part of on-line adaptive control 
or media replacement strategy for AFM, once they have been correlated with existing 
abrasive wear signals. 

AFM process is very much in its early stages of development. Sustained research is 
required to transform the process into mature technology, for its successful induction into 
the industiy. This thesis makes an attempt in this direction. The scope of the present work 
and organiaation of the thesis are described in the following sections. 



16 


Introduction and Literature Review 


1.5 Scope and Objectives of the Present Work 

Although the investigations carried out by various researchers are very good; they 
seem to be rather lacking in theoretical treatment. Both theoretical and empirical studies of 
abrasive flow machining are greatly hampered by the inherent random nature and 
multiplicity of variables. Not much information is available in the literature, which deals 
with theoretical analysis of material removal mechanism in AFM. The relationship 
between process parameters and performance characteristics of the process are not 
completely known. The exact mechanism by which the individual abrasive particles 
accomplish material removal is only partially understood. There is a need to carry out 
systematic theoretical analysis of viscous media flow to determine the stresses and 
machining forces acting on the workpiece. This will help in modeling and proper 
understanding of the mechanics of the machining process during AFM. 

The following are the specific objectives of the present work: 

• To conduct experiments to study the effects of various process parameters on material 
removal and surface roughness. 

• To analyze the flow of media in AFM process using finite element method by 
considering the media as viscoelastic material. 

• To develop a model to determine the material removal and surface finish in AFM 
process and to validate it by comparing with the data obtained by experimentation. 

• To develop a model for the topography of media and estimate active grain density by 
microscopic technique, and to simulate the random surface generated in AFM process. 

• To analyze the specific energy and thermal aspects of AFM process. 

• Modeling of AFM process by neural network, and optimization using neural network 
and genetic algorithm. 


1.6 Organization of the Thesis 

The total thesis write-up is divided into seven chapters and four appendices. First 
chapter of the thesis comprises introduction of traditional finishing methods and their 
limitations, introduction and importance of abrasive flow machining process, and review 
of literature. Based on the literature review, the scope and objectives of the present work 
are presented. 
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Description of the experimental set-up, details of experimentation and the 
associated results constitute chapter 2 of the thesis. Experiments have been conducted on 
the AFM machine designed and fabricated at Indian Institute of Technology, Kanpur. 
Experiments have been planned using central composite rotatable design to get the useful 
inferences by performing minimum number of experiments. A parametric study of process 
performance is also carried out. 

Chapter 3 presents a model for analysis of media flow. The method involves the 
use of finite element technique to solve the governing equations in terms of primary 
variables (velocities and pressure). The theoretical analysis of media flow predicts the 
stresses developed during the flow. 

Modeling for material removal and surface roughness using abrasion theory is 
also described in chapter 3. The model is based on accumulated plastic flow of material by 
repeated indentation by moving abrasive particles. The model includes the operating 
conditions of the process e.g. abrasive size, extrusion pressure, workpiece material 
hardness and size of the workpiece to be machined. 

A stochastic simulation method using random sampling technique is proposed to 
evaluate the active grain density in chapter 4. A microscopic technique has also been used 
to estimate the number of active grain density at various mesh sizes and concentrations. 
The modified model to simulate the random surface generated and material removal in 
AFM is also discussed. 

The chapter 5 is devoted to forces, specific energy and thermal aspects of AFM. A 
model has been proposed for the determination of specific energy and tangential forces 
in AFM process. A simple model of heat transfer in AFM has also been developed that 
considers heat flow to the workpiece and abrasive media. 

Chapter 6 deals with the modeling of AFM by neural networks. The thesis presents 
an architecture of back-propagation neural network model for prediction of material 

removal rate and surface roughness in AFM. 

The neural network has been used for optimum selection of input parameters of 
AFM process. Optimization of the AFM process has also been performed using genetic 
algorithm (GA). The optimization results of neural network are compared with the results 

obtained by genetic algorithms. 
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The conclusions drawn from the present work are enumerated in chapter 7. 
Possible avenues for future research in AFM are also outlined. The references and 
appendices are given at the end of the thesis. 


Chapter 2 
Experimentation 


This chapter presents the effects of AFM process parameters on process 
performance determined experimentally. Experiments have been conducted on mild steel 
workpieces to determine the material removal and change in surface roughness under the 
specific machining conditions. The experiments have been planned according to central 
composite rotatable design. The effects of percentage concentration of abrasives, extrusion 
pressure, reduction ratio and number of cycles on material removal and surface roughness 
have been discussed. The following sections describe the experimental set-up and details 
of experimentation. 

2.1 Working of AFM Set-up 

The Abrasive Flow Machining set-up has been designed and fabricated at 
Manufacturing Science Lab., Indian Institute of Technology Kanpur, keeping in view the 
fundamental mechanism of the process and the basic requirements for research purposes. 
Media composed of silicon carbide abrasives (mesh size = 50-60) and polyborosiloxane 
carrier was used (Supplied by Extrude Hone Limited, U.S.A). Polyborosilixane carrier is 
soluble in acetone. Abrasive grains can be separated from media by dissolving the carrier 
in acetone. In AFM, media is extruded to and fro through the passage formed around the 
workpiece by utilizing two opposed piston media cylinders. The passage around the 
workpiece surface to be finished is formed by means of fixture arrangement. 

The Abrasive Flow Machining set-up shown in Figures 2. 1 and 2.2 consists of the 
following main components [24] 

1 . Media cylinders 

2. Pistons 

3. Hydraulic drive and controls 

4. Workpiece fixture 

5. Frame and housing. 
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Media cylinders 

The main function of the media cylinder (Fig. 2.2) is to embody the media at 
required extrusion pressure. It should also guide the piston movement. It can withstand 
pressures upto 100 bar. As the media is moving back and forth, the abrasion will take 
place on the cylinder inner surface. Hence, replaceable cylinder liners are also provided. 

Pistons 

The piston acts to transmit the extrusion force applied by the piston rod of 
hydraulic cylinder to the abrasive media. It is used to extrude the media through the media 
cylinder to the workpiece. The main extruding surface is formed by the piston head, which 
accommodates two cast iron rings to provide sealing and reduce frictional resistance. 
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Hydraulic drive and controls 

In the AFM, the operations are powered by the hydraulic drive, 'i’hc hydraulic system 
refers to the complete assembly of component parts that transmit and control the fluid 
power. Fig. 2.3 shows the block diagram of a hydraulic system, big 4 shows the 
hydraulic circuit which has following components. 

1 . Reservoir (R) 

2 Strainer (S) 

3 Variable delivery pump (VDP) 

4. Flexible coupling 

5 Electric motor (M) 

6 Connectors 

7 Relief valve (RV) 

8. Check valve (CV) 

9. 3-position 4-ways direction control valve (DCV) 

10. Single acting cushioned cylinders (HCL & HC2) 

The pump VDP withdraws oil from the inlet line C and delivers it to the outlet line 
D at a volumetrically uniform rate for any particular setting. The direction of flow is 
reversed by means of valve DCV so that line P is connected to either of the cylinder 

To ensure that the maximum pressure is limited, a shunt line G is provided, 
communicating with reservoir R by way of relieve valve RV3. It is apparent that, with this 
arrangement, the rate of advance of the piston Pi is determined by the rate of withdrawal 
of oil from the hydraulic cylinder HC2. As oil is passed through the pump VDP at a 
volumetrically uniform rate, the piston would advance at uniform speed provided that the 
oil is incompressible. 

During forward stroke of the cycle, oil enters from the port P of the direction 
control valve (DCV) and flows into the line A connected to check valve and relief valve 
for pressure control in the cylinders. Initially RVl and RV2 are adjusted manually for a 
particular pressure. In forward stroke the oil passes from CVl into HCl till the pressure in 
the hydraulic cylinder HC2 reaches the adjusted pressure. As the pressure exceeds the 
adjusted value, the RV2 gets opened and oil comes out from HC2, making piston PI 
traversing in downward direction. After completing forward stroke, the direction of DCV 
is reversed by relays and oil flow is reversed to complete one cycle. 
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Workpiece fixture 

The function of fixture or tooling is simply to hold the workpiece in position and to 
direct the flow of media to pass around the surface to be finished The design of the 
workpiece fixture is different for different shapes of workpieces. Keeping in the view the 
fixture design considerations, fixturing arrangement, as shown in Fig. 2.5 is fabricated to 
fulfill the requirements The fixture is fabricated for finishing internal surface of 
cylindrical workpiece The workpiece fixture consists of main body, which 
accommodates the workpiece and tooling The tooling is different for different internal 
diameter of workpiece to direct the abrasive media at a particular angle 



tooling 

workpiece 

hole for themocouple 


Fig. 2.5 Schematic diagram of fixture used for AFM 
Frame and housing 

Frame and housing are used to accommodate media cylinders, hydraulic cylinders, 
workpiece fixture, etc together as a single unit. It is capable to withstand all kind offerees 
acting on it. 

2.2 Procedure and Precautions to Conduct Experiments 

The following procedure should be followed during experimentation on AFM set-up. 

1 Abrasive media is filled tightly in media cylinder. It is preferred to fill the media in 

t 

bottom cylinder, keeping piston in extreme bottom position. 
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2 Initial weight and surface roughness of the workpiece are measured before machining 

3 Then, the workpiece is fixed in the fixture as shown in Fig 2 5 and the media is filled 
in the workpiece fixture 

4 The top media cylinder assembly is lifted by the screw arrangement and the fixture 
(filled with media) is kept over the bottom media cylinder as shown in Fig 2 2 The 
fixture is tightened, so that no media should leak out during experimentation, 

5 The motor is started, ensuring that the pressure is at 0 bar 

6 The voltage is adjusted to approximately 6V in the DC power supply and then, the 
control mechanism for solenoids should be switch on to change the direction of flow 
of oil. 

7. Keeping the control switch in neutral position, the pressure is increased to the required 
pressure by the pressure valve arrangement. 

8 After this, the automatic control switch for solenoids is kept in on position and it 
automatically changes the direction of flow of the oil. One AFM cycle is completed 
when the media is extruded from the lower cylinder to the upper cylinder and back 
again to the lower cylinder. 

9 After completing required number of cycles, the automatic control switch and motor 
are switched off. 

10. After machining, the workpiece fixture and workpiece are taken out from AFM set-up, 
and the final weight and surface roughness of workpiece are measured 

The following precautions should be observed while conducting experiments: 

1 . Before starting the experiments, it should be made sure that the control switches are 
working properly. 

2 Fixture should be placed in alignment to the media cylinder. 

3 . The plates should be hold very tightly above the fixture, so that no media leaks out 

4. Both manual and automatic switches provided for control of solenoids should not be 
kept ‘on’ at the same time. Otherwise it may damage total solenoid circuit. 

5. Design of workpiece fixture should be such that the workpiece is held tightly. The 
workpiece should be cleaned after machining, but before weighing it. 
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2.3 Guidelines to Design Fixture / Tooling 


The basic functions of the fixture / tooling for AFM process are as follows 

• It should assure high accuracy of the parts. 

• It should position the workpiece at the desired place and hold it tightly. 

• It should guide media, and form the restriction passage so that the finishing takes place 
in the desired manner 

• If possible, it should provide interchangeability. 

To fulfill the above basic functions, fixture should contain the following components / 
elements: 

• A sufficiently rigid body to withstand all types of loads coming on it. 

• Locating elements to position the workpiece at a required place. 

• Clamping elements to hold the work tightly. 

• Elements for positioning a fixture on the AFM set-up. 


2.4 Design of Experiments 

Experiments conducted were planned using statistical techniques so that we get 
useful inferences by performing minimum number of experiments. If the yield or response 
‘y’ is a function of the levels of quantitative variables then we may write [25], 

yn = ^ (Xiu, X2u, , Xku) + eu (2. 1) 

where, u= 1,2, ,N represent the N observations in the factorial experiment, and Xiu 
represents the level of the ;th vector in the wth observation 

The function (j) is called the response surface. The residual eu measures the 
experimental error of the u* observation. A knowledge of the function ^ gives a complete 

summary of the results of the experiments and also enables us to predict the response for 

the values of the x^ that are not tested in the experiment. When the mathematical form of 

1 

<f) is not known, this function can sometimes be approximated satisfactorily, within the 

experimental region by a polynomial equation in the variables Xiu. The important criteria in 
selecting design for fitting a second order response surface is that the computations should 
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not be too difficult Considering desirable properties, Box and Hunter proposed the criteria 
of rotatability. In the experiments conducted for AFM process analysis, the Central 
Composite Rotatable Design [25] is used The general from of the quadratic polynomial 
fitted to the experimental data is given below 


k k k 

yu = bo+ Eb,x,„+ Zb„x,;+ ZbyX.^x 

1=1 1=1 Kj 


ju 


( 2 , 2 ) 


where, k is no of variables and u is number of observation 


From the results of any experiment, we can compute the standard error of yu, at any 
point on the fitted surface The standard error will be a ftinction of the co-ordinate x,u of 
the point In a rotatable design, this standard error is the same for all points that are at the 
same distance from the center of region This property is reasonable one to adopt for 
exploratory work, in which the experimenter does not know in advance how the response 
surface will orient itself with respect to the x-axis Consequently he has no rational basis 
for specifying that the standard error of yu should be smaller in some directions than the 
others 

The design for four x - variables (say, xi= extrusion pressure, X2 = number of 
cycles, X3 = reduction ratio, and X4 = percentage concentration of abrasives) is shown in 
the appendix A (Table A.1). The columns headed xi, X2, X3 and X4, which specify the 
actual combination to be used, constitute the plan of experiment. 

Subsequent steps are as follows. 

1 . The columns headed xo, xi, X2, X3, X4, xi^, X2^, xj^, x/, X1X2. X1X3, X1X4, X2X3, X2X4, X3X4 
as shown in Table A.l (Appendix A) are completed. The two way array with 16 
columns and 3 1 rows comprises the X - matrix of the x-variables The corresponding 
values of the response y are placed on the right. 

2. The sum of products of each column in the matrix X with the column of y values are 
calculated. These sum of products are denoted by ( 0>' ), (!>') , ( 2 y) and so on. 

3 From the values of (O3; ), (ly) , etc the regression coefficients are computed directly by 
the equations given below [25], The only auxiliary quantity needed is YjO^y) > found 
by adding the crossproducts of all the squared terms with y . 



28 


Experimentation 


Z>o=0 142857(0;;) - 0 035714X«> 

Y^iiy - (1 \y) + (22>') + (33;;) + (44;;) 

b =0 041667 (ly) (2.3) 

b„ =0.03 1250(/7;;) + 0.003720£(//y) - 0 035714(0;;) 
b,j=0 0625(//» 



Fig. 2.6 Dlustration of AFM process 

In the present analysis of the experiments, the effect of four variables viz. pressure 
Pg , number of cycles N , percentage abrasive concentration by weight C , and reduction 

ratio i?gare studied on material removal and change in surface finish. Reduction ratio 

[ = 1 - (r^ /r/ ) ] expresses the reduction in cross sectional area of flow of media. It is 

defined as the difference between cross sectional area of media cylinder and workpiece, 
divided by the cross sectional area of media cylinder. Here, is internal radius of media 

cylinder and is inner radius of cylindrical workpiece (Fig. 2.6). Slug length (this 
term is used in subsequent sections) varies with reduction ratio. It is defined as the length 
of media passed through workpiece in one stroke ([/„, = /^ / (1 - 7?^)] . Here, is length of 

stroke. Percentage concentration of abrasives by weight C can be defined by the 
following equation. 

^ ^ weight of abrasives 

V total weight of abrasives and carrier j 


*100 


(2.4) 
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A preliminary step is to set-up the relationship between the coded x - scales and the 
original scales in which the levels are recorded In design scale, the lowest and highest 
values of x are -2 0 and 2.0 So we take, 

X = -2 0 when = 20 bar; N =5 cycles, C = 56 %, R^=0U 
X =2 0 when = 80 bar, N =25 cycles, C = 76%, = 0 97 

Then, 

X = a + b * (variable) 

'a' and 'b' are chosen to satisfy the desired conditions at the end of the scale 

After substituting the lower and higher values of variables at original scale 
corresponding to -2 0 and 2 . 0 , the following expressions, which describe the relation 
between x and original scale of variables, are obtained' 

Xi = -3.3333 +0.0666 p^ 

X 2 =-3.0+02 A 

X 3 = -27.8461538 + 0 3076923 R, 

X 4 = -13.2 + 0.2 C 

where xi, X 2 , X 3 and X 4 denotes the values of x corresponding to pressure, number of 
cycles, reduction ratio and percentage concentration of abrasives respectively. 

From these equations, extrusion pressure, number of cycles, reduction ratio and 
percentage concentration of abrasives in media corresponding to level x = - 1 , 0 , 1 are 
determined. 

The complete conversation is tabulated in Table 2.1. 


Table 2. 1 Conversion Table 


X 

Extrusion 

pressure, 

bar 

No. of 

cycles 

Reduction 

Ratio, % 

Concentration 

% 

- 2.0 

20 

5 

84.00 

56 0 

- 1.0 

35 

10 

87.25 

61.0 

00 

50 

15 

90.50 

66.0 

1.0 

65 

20 

93 75 

71.0 

2.0 

80 

25 

97 00 

76.0 
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With the above conversions, the experiments have been conducted according to the 
order given in the Table 2.2. The table also gives the responses - material removal and 
change in surface roughness observed during experimentation. The experiments 

with lowest abrasive concentration have been conducted first, then in the increasing order 
of percentage concentration of abrasives, due to limited availability of media. For 
increasing percentage concentration of media, abrasives were mixed in the same media 
used in the previous experiments. Specifications of workpiece, abrasives and carrier used 
during experimentation are given in Appendix C. 


Table 2.2 Plan of experiments 




Designed 

parameters 


Responses 

Ex. 

Pressure, 

Number of 

Reduction 

Concentration 

Material 

AT?, 

No. 

bar 

cycles 

ratio % 

% 

removal, g 

(itm) 

B 

50.0 

15 

90.50 

56 

0.0147 

0.43 

2 

35.0 

10 

87.25 

61 

0.0090 

0.35 

3 

65.0 

10 

87.25 

61 

0.0146 

0.50 

B 

35.0 

20 

87.25 

61 

0.0177 

0.52 

5 

65.0 

20 

87.25 

61 

0.0241 

0.77 

6. 

35.0 

10 

93.75 

61 

0.0123 

0.56 

7. 

65.0 

10 

93.75 

61 

0.0167 

0.69 

8. 

35.0 

20 

93.75 

61 

0.0214 

0.73 

9. 

65.0 

20 

93.75 

61 

0.0297 

0.94 

10. 

20.0 

15 

90.50 

66 

0.0075 

0.32 

11. 

80.0 

15 

90.50 

66 

0.0282 

0.74 

12. 

50.0 

5 

90.50 

66 

0.0065 

0.31 

13. 

50.0 

25 

90.50 

66 

0.0274 

0.77 

14. 

50.0 

15 

84.00 

66 

0.0063 

0.36 

15. 

50.0 

15 

97.00 

66 

0.0495 

1.01 

16. 

50.0 

15 

90.50 

66 

0.0217 

0.61 

17. 

50.0 

15 

90.50 

66 

0.0229 

0.54 
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18 


15 

90.50 

66 



19. 

50.0 

15 

90 50 

66 

0.0211 

0.59 

20. 

50 0 

15 


66 



21 

50.0 

15 

90,50 

66 



22 

50 0 

15 

90.50 

66 



23, 

35 0 

10 

87.25 

71 



24 

65.0 

10 

87 25 

71 

0 0186 

0 57 

25 

35 0 

20 

87.25 

71 

0 0197 

0 65 

26 

65.0 

20 

87 25 

71 

0.0285 

0.89 

27 

35 0 

10 

93.75 

71 

0.0209 

0.54 

28. 

65 0 

10 

93.75 

71 

0.0335 

0.71 

29 

35.0 

20 

93.75 

71 

0.0363 

0.79 


65.0 

20 

93.75 

71 

0.0465 

1.03 

31 

50.0 

15.0 

90.50 

76 

0 0279 

0 85 


2.5 Material Removal Measurement 

The weight of the workpiece was measured using the available electronic weight 
balance “AFCOSET FX-400” having accuracy of O.OOOlg. It is essential to clean the 
workpiece properly before measuring the weight By measuring the weight of the samples 
before and after machining, the material removed has been calculated by taking their 
difference. The data of material removal measurements are given in Table A.2 (Appendix 
A) and their average values in Table 2.2. 

2.6 Surface Roughness Measurement 

To study the change in surface texture due to various process parameters, the 
surface roughness measurements are essential. For this purpose, the instrument used was 
"TAYLOR HOBSON SURTRONIC 3P " having least count of 0.01pm. It gives on its 
digital scale the , Ry and Rf„ values. The improvement in surface finish has been 

determined by measuring R^ , Ry and R^^ values before and after machining. 
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is most widely used description of a surface. For determining value over a 

sampling length, a mean line is computed such as the area of the solid above the line is 
equal to the area of the void below. The arithmetic average of the deviation up and down 
from this theoretical mean line is R^ 

Ry is the distance from the highest peak to the deepest valley Ry is useful for 

detecting the finishing problems, like pick-up or areas which have not cleaned up 

R^^ is also called Rz DIN. For determining R ^^ . the measured length is divided 
into five equal lengths and then average of the highest peak to the deepest valley of each 
of the five lengths is equal to Rj„ It is likely to be a slightly smaller than Ry because one 

deeper is diminished by the finer finish of the other four values 

The data of surface roughness measurements of workpiece surface before and after 
machining are given in Table A 3 (Appendix A). The data representing improvement in 
surface roughness are given in Table A.4 (Appendix A). The average of change in vSurface 
roughness (A/?^) for each experiment is given in Table 2.2. 


2.7 Response Surface Analysis (RSA) 

The equation of general quadratic response surface for four variables is of the form 
y = bo + bixi + b2X2 + bsxa + b4 X4 + bnXi^ + b22 X2^ + bjs -i- b44 X4^ + bi2 X1X2 

+ bi3 X1X3 + bi4 X1X4 + b 23 X2X3 + b 24 X2X4 + b34 X3X4 (2.5) 

where y is response under study, bo, bi. .... , b34 are constants, and xi, X2, xa and X4 arc 
process variables as mentioned in section 2.4, The quadratic response surface for each 
response is obtained by calculating constants bo, bi, , b34 of equation ( 2 . 5 ). 

2.7.1 Material removal 

The constants bo, bi, , b34 are calculated from eqn. 2.3 using the responses of 

table A.2 (Appendix A) or table 2.2. The following response surface equation is obtained 
for material removal (MR). 

MR = 21.714 + 4.358 Xi + 5.350 X2 + 6.658 X3 + 4.025 X4 - 0.703 Xi^ - 0.928 X2^ 

+ 1 808 X3 + 0.158 X4^ + 0.262 X1X2 + 0.487 X1X3 + 0.862 X1X4 + 0.900 X2X3 

+ 0.375 X2X4 + 2.750 X3X4 (2.6) 
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2.7.2 Surface roughness 

The constants bo, bi, , b 34 are calculated from eqn 2.3 using the responses of 
change in surface roughness given in table A 4 (Appendix A) or table 2 2. The following 
response surface equation is obtained to predict a change in surface roughness (Ai?^) 
generated by AFM process 


AR^ =0.617 + 0 lOlxi + 0 119x2 + 0.107x3 + 0 067x4-0.009x1^-0.006x2^ 

+ 0 029 X 3 ^ + 0.006 X 4 ^ + 0.023 xiX 2 - 0 0006 xiX 3 + 0.0018 X 1 X 4 + 0.0068 X 2 X 3 
+ 0 0093 X 2 X 4 - 0.01 18 x 3 X 4 


(2.7) 


2.8 Results and Discussion 

Responses have been calculated using response surface equations (2.6) and ^2 7), 
which have been obtained from experiments based on central composite rotatable design, 
to study the effects of various parameters on material removal and change in surface 
roughness. In the following chapters, the results obtained using these eqns (2.6) and (2.7) 
are termed as response surface analysis (RSA) results. Figures 2.7 and 2 8 show that the 
experimental points are close to those predicted by RSA. Hence, RSA results can be 
considered as representatives of experimental results At the center point of the 
experimental design, eight experiments were carried out for same value of variables. As 
there is very little difference in the results for material removal and surface roughness 
values obtained during these eight experiments (Table 2.2), one can conclude that process 
repeatability is good. 

2.8.1 Material removal 

Figure 2 7 shows the effect of percentage concentration of abrasives by weight on 
material removal. As the percentage concentration of abrasives in media increases, 
material removal also increases, because higher concentration of abrasives results in 
higher number of active grains per unit area and hence higher material removal. Further, 
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grains can sustain larger cutting forces due to increase in viscosity of media at higher 
concentration. 

It is found that material removal increases with increase in extrusion pressure (Fig 
2.8). Depth of indentation increases with increase in extrusion pressure, which results in 
higher area of indentation of a grain and hence higher material removal occurs 

Material removal also increases with increase in reduction ratio (Fig. 2.9), or due 
to increase in slug length of media passed per unit time. As volume of media passed across 
the workpiece in an AFM stroke is fixed, slug length of media passed across workpiece 
becomes more at low value of cross sectional area of workpiece and vice versa. At higher 
reduction ratio, the pressure acting on the workpiece surface is also higher, leading to 
higher depth of indentation. Both these factors result in higher material removal at higher 
reduction ratio. 

Material removal value increases non-linearly with the number of cycles (Fig 2.10). 
The rate of material removal slightly decreases with the increase in number of cycles. This 
can be explained as follows Initially unfinished surface has sharp peaks, the peaks get 
machined and become comparatively flater after a few cycles of machining. 



Fig. 2.7 Variation of material removal with percentage concentration of abrasives 

(/?,= 50 bar, i?, = 0.905, N =15 cycles) 
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Fig. 2.8 Variation of material removal with extrusion pressure 
{R^ = 0.905, N =15 cycles, C = 66%) 



Fig. 2.9 Variation of material removal with reduction ratio 
(/ 7 g= 50 bar, N =15 cycles, C = 66%) 
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Fig. 2.10 Variation of material removal with number of cycles 
(/?, = 50 bar, R, = 0.905, C = 66%) 

2.8.2 Surface roughness 

The effect of percentage concentration of abrasives on the change in surface 
roughness indicates the increase in change in surface roughness with increase in 
percentage concentration of abrasives (Fig. 2.11'). With higher concentration of abrasives, 
the number of abrasives taking part in machining will be more^ resulting in more abrasion, 
and hence surface finish obtained for a specified number of cycles will be better. It is also 
observed that with increase in percentage concentration of abrasives, the improvement in 
surface roughness value is more rapid. 

Figure 2.12 exhibits the role of extrusion pressure on surface roughness value 
obtained after the specified number of machining cycles It is observed that change in 
surface roughness value increases with increase in extrusion pressure but its rate 
decreases. 

Figure 2.13 shows the effect of number of cycles on change in surface roughness 
value of the workpiece surface. As the number of cycles increases, the roughness value 
decreases. In other words, change in surface roughness value increases with number of 
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cycles (Fig 2 13) Initially, rate of improvement in surface roughness value is higher 
because of availability of peaks in early stages, but as these peaks get machined, the rate 
of improvement in surface roughness value slowly decreases 

Figure 2 14 illustrates the effect of reduction ratio on change in surface roughness 
value of workpiece surface It is observed that change in surface roughness value increases 
with increase in reduction ratio, at specified number of cycles At higher reduction ratio, 
the material removal is higher due to higher slug length of flow and higher normal forces 
As a consequence, peaks get abraded more in specified number of cycles, and hence, 
surface roughness value is lowered 



Fig. 2.11 Variation of change in surface roughness with percentage 
concentration of abrasives 
(/ 7 g= 50 bar, = 0.905, N =15 cycles) 



Change In surface roughness, Change in surface roughness, 

fdm H fxm 
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{R^ = 0.905, N =15 cycles, C = 66%) 



5 10 15 20 25 

Number of cycles 


Fig. 2.13 Variation of change in surface roughness with number of cycles 

50 bar, = 0.905, C = 66%) 
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Fig. 2.14 Variation of change in surface roughness with reduction ratio 
(jPg= 50 bar, N =15 cycles, C = 66%) 

2.9 Conclusions 

AFM experiments on mild steel workpieces using SiC abrasives mixed with 
polyborosiloxane carrier indicate that required surface finish of workpiece surface could 
be achieved by controlling the AFM process parameters. Central composite rotatable 
design made possible to study the effect of four important process parameters on process 
performance by performing only 31 experiments. A comparison between RSA results 
(derived from equations (2 6) and (2.7)) and experimental data reveals that RSA results 
can be considered as representative of experimental results. An increase in extrusion 
pressure increases the depth of indentation of abrasive and increase in percentage 
concentration of abrasives increases active grain density and media viscosity, hence 
material removal and change in surface roughness are higher at higher extrusion pressure 
and concentration of abrasives in media. Material removal and change in surface 
roughness also increase with increase in reduction ratio and number of cycles. 




Chapter 3 

Modeling of Material Removal and Surface Roughness 
in AFM 


3.1- Introduction 

Not much information is available in the literature, which deals with the 
theoretical analysis of material removal mechanism in AFM. The relationship between 
process parameters and performance characteristics of the process is not known 
completely. The exact mechanism by which the individual abrasive particle 
accomplishes material removal is only partially understood. There is a need to carry out 
systematic theoretical analysis of the media flow to determine the stresses and 
machining forces acting on the workpiece. This will help in the modeling and proper 
understanding of the machining action during AFM. Once machining forces acting on 
the workpiece surface, developed due to media flow, are known, an expression can be 
established for material removal and surface roughness in terms of extrusion pressure 
and other process variables. The numerical simulation of viscous flow has been the 
topic of many researchers [26 -30]. Early numerical formulation of viscous materials 
involved finite element technique based on the Galerkin models, and finite difference 
technique based on the stream-function vorticity formulations. Zienkiewicz et al. [26] 
have presented the flow-formulation approach in forming and extrusion, investigating 
two techniques, viz. the pressure-velocity formulation with Lagrangian constraints and 
the penalty-function approach. Dixit et al. [27] carried out deformation analysis of the 
steady-state wire drawing process and Reddy et al. [28] carried out analysis of tube 
extrusion by FEM using the mixed pressure-velocity formulation. The initial 
developments of finite difference and finite element methods used for solving non- 
Newtonian flow have been reviewed by Crochet et al. [29]. Rajeshwar et al. [23] 
presented a simulation model to determine the characteristics of media flow during 
machining. The finite difference method was chosen for obtaining the solution. Lim 
and Dunne [31] have performed finite element modeling of the axisymmetric extrusion 
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of an aluminum matrix composite containing silicon carbide particles using continuum 
level constitutive equations In their analysis for making the problem mathematically 
tractable, they considered metal matrix composite material containing high volume 
fraction of reinforcement as isotropic and homogeneous. Wilson et al [32] reported that 
if the settling velocity of discrete solid particles transported by a carrier liquid is very 
small, then the slurry can be considered as single phase. If the mixture can be 
transported with a uniform solid concentration across the pipe, then during such 
transport the mixture can be viewed as homogeneous. The mixture behaves as 
Newtonian at low solid concentrations and non-Newtonian at higher concentrations. 

In the present case of AFM, the physical state of the media is that of a very 
viscous fluid, and able to adopt the shape of any container. The state of abrasives in the 
carrier is very similar to what has been stated in the preceeding paragraph through 
diversified examples. Hence, for simplification of the analysis, it is assumed that the 
abrasives are uniformly distributed in the carrier (polyborosiloxane), so that the media 
acts as a homogeneous mixture (continuous phase) [20]. Others [23] also have analysed 
the deformation of media in AFM process considering it as homogeneous and isotropic. 

The objective of the present chapter is to establish a model for media flow in 
AFM, and hence to determine the stresses developed during the flow. These computed 
stresses can be used to evaluate the machining force and hence, material removal and 
surface roughness can be obtained in AFM for the given machining conditions. This 
chapter also describes the simulation of flow of media as a viscoelastic material in 
AFM process. Solving this problem using available classical methods would be very 
difficult and tedious. It is also not possible to apply the classical methods to complex 
shapes like turbine blades which are finished by AFM process. Hence, finite element 
method has been used to solve such complex geometry problems. Finite element 
technique is used to solve the governing equations in terms of primary variables 
(velocity and pressure). Further, theoretical approach to the estimation of material 
removal and surface finish during abrasive flow machining process has also been 
presented. 
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3.2 Analysis of Media Flow in AFM Process 

The mathematical representation of the flow of media in APM process involves the 
equations of continuity and momentum, and the constitutive equations The following 
assumptions are made to simplify the analysis: 

1 The media characteristics are similar to a Maxwell fluid [33] 

2. The media is isotropic and homogeneous. The media properties are independent of 
temperature, and constant with time and space. 

3 For cylindrical workpiece, media flow is taken to be axi-symmetric. 

4. The media flow is assumed to be steady and fully developed 

Steady state, incompressible and axi-symmetric form of the continuity, momentum 
and constitutive equations is given below. 


3.2.1. Basic equations 

(A) Continuity equation (conservation of mass) 

For steady state, incompressible and axi-symmetric flow of media, neglecting the 
volume of material removed* from work surface, the condition of volume constancy can 
be expressed as, 

Err + see +£'22 =0 (31) 


where, £rr, see and e-a. are the components of strain rate tensor in r, (9, and z 
directions, respectively and can be expressed in terms of the velocity Vj as 


+ ^ ( 3 .^) 

2[<^Xy, 

(B) Momentum equations (balance of momentum) 

The momentum equations for a steady-state axisymmetric flow of media 

(neglecting the body forces) have the following form: 

(3.3) 


r j. 

P [ V, ^ + V 
dr 


. dp p 1 dip T^) 


+ :£Es-] = o 


dz 


and 


• The total volume of material removed during AFM is extremely small, because it is a finishing process and 
not a bulk material removal process 
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= 0 (3.4) 

dr dz dz r dr dz 

where, p is density of media, p is pressure (hydrostatic part of stress tensor), and is 
the extra-stress (deviatoric part of stress tensor) 

(C ) Constitutive equations 

The media exhibits both viscous and elastic properties [33], In the present analysis, 
simple constitutive equation which contains a minimum number of material parameters 
such as model of Maxwell is selected According to Maxwell’s model, the constitutive 
equation is given as follows- 


+ X 

dt 


where, X is the relaxation time and p is viscosity coefficient. 


— is material time derivative For steady state viscoelastic flow, the material derivative 
dt 

of the stress deviator is defined as 


where v,^ are components of velocity tensor. 

The expanded form of the eqn (3 5) is given by the following equations. 




+ 2 (Vr ^ + Vz -^ ) = p 


r„+ 1 (Vr ^ + Vz ^) = p 

dr d z 


^ (Vr -T^+ Vz -^)= p s. 


(3.10) 
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For an incompressible fluid, the total stress tensor a ^ is given by 
P 5, + (3 11) 

where, are components of the identity tensor ( or Kronecker delta) 


3.2.2 Non-dimensionalization 

Non-dimensionalization is usually required to avoid ill-conditioning and numerical 
difficulties. The non-dimensionalization of all physical quantities used above has been 
done using characteristic dimensional quantities as given below 



and W = A ^ (3 12) 

r^, 


where, W is non-dimensional parameter which is sometimes called Weissenberg 


number, and ~ - -7- ; is the flow stress of material; A-,„is the radius of 


V /r > 


workpiece, and Vp is velocity of media flow at inlet 


Substituting non-dimensional parameters in continuity, momentum and 
constitutive equations (3.1, 3.3, 3 4 and 3.5), we obtain 

0 

^ ^] = 0 (3.14) 

Po ^ ^ d r r d r r d z 




(3.15) 
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dr 

r„ + W ^ 




(3.16) 


dt ' 

For convenience in writing, the overbar (-) from non-dimensional parameters is 


dropped in the succeeding description. Since Reynolds number 




is found to be 


Po 


typically of the order of 10'* for all the cases considered here, the inertial terms of eqns 
(3 14 and 3.15) are ignored in the next section. 


3.2.3 Galerkin or weak formulation 


The finite element formulation of flow of the rigid plastic material in terms of 
primary variables (pressure and velocity) has been presented by Dixit et al. [27]. Similar 
formulations with some modifications are possible for the flow of AFM media, which is a 
viscoelastic incompressible material 


In the present analysis, the Galerkin weighted residual method [34] is employed to 
reduce the governing differential equations by using appropriate weighting functions for 
the momentum and continuity equations. Let Vr, Vz, and p constitute the functions that 
satisfy all the essential boundary conditions. Then Vr, Vz, and p constitute a weak solution 
if the following integral equation is satisfied. 


J [ { ^rr + See + ffzz } M/p 
a 


+ { 


^ P 
d r 


\ ^rr) + lee. 
r dr r 


dr, 


d z 




+ { ££. 1 

d z r dr 


— — } ] 2 n r dr dz 

d z 


(3.17) 


where. M/p, and m/^ are the weight functions which satisfy the homogeneous versions 
of the boundary conditions and Q represents the area of the domain. 

Performing the integration by parts on second and third segments of the above 
equation, the following weak form is obtained: 


J [h +I 2 ] 2 ;>r r dr dz - j I3 dTj - j I4 oTj = 0 
^ r, fj 


(3.18) 
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where, Ii = + £ 2 z}m> 


(3 19) 


h--p Seei^)^ ]+ + + T^{^ + ?^)} 

oz or r or oz 


I 3 = t, 


(3 20) 
(3 21) 


I4 = t,. 


(3 22) 


and r, and Fj are the boundaries where the traction components t,. and are specified 
The term (w), Sgg (w ), e^{w) and £^(m>) are defined as 


<5* 


d w. 


or a z 




, £ee{'^)= — \ and£,^(w)=- 

r 2 


1 w,. , d-w. 


dr d z j 
(3 23) 


3.2.4 Finite element approximations 

In the present analysis, Langrangian biquadratic shape function for Vr and Vz , and 
bilinear shape function for p over a 9 noded quadrilateral element are used The 
approximations for Vr and Vz are as follows : 

Kr 


iV, 0 iVj 0 A^9 0 

0 A^i 0 A^2 0 iVg 


= [iVv]{vn 


(3 24) 


(v^r 

(vz,r 

In Galerkin formulation, the weight functions corresponding to velocities are 
approximated using the same shape functions as that for velocity. Therefore, 

{“'-}= {“' I “ 


(3 25) 
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The vectors {v®} and (w®} contain the velocities at the nodes of elements and their 
corresponding weight functions respectively. The approximation for p is as follows: 


p = [n( Ni Ni N(] 


Pi 

pI 

p! 

pI} 


(3.26) 


i.e. p= {N^fip^ } 


(3.27) 


The weight function corresponding to pressure is approximated using the shape functions 
as that for pressure. 

= {If/ {«.;) (3.28) 

The deviatoric stresses are approximated by the same shape functions as used for pressure. 
Therefore, 

[r] = [NMr^] (3.29) 


0 




0 

0 0 
0 0 


0 Q 0 Q 0 . 
0 Q 0 0 0 . 

‘JlNy 0 0 0 V2A^2 0 • 

0 iVj 0 0 0 . 


iV4 0 0 0 

0 0 0 

0 0 V1/V4 0 
0 0 Q 


(3.30) 


'(V)i 

(^zz)l 



(^rr)4 

(^2r)4 

(^^2)4 


(3.31) 



Modeling of Material Removal and Surface Roughness 


49 


3.2.5 Finite element equations 

To facilitate derivation of the finite element equations, it is first necessary to 
express eqn (3 1 8) in vector form The following vectors are defined for that purpose 
The strain vector is given by 


^60 


dr 
d z 


1 


4i 

V. 


dv^ d\. 


dr d z 


(3 32) 


Substitution of the approximating expressions in the above equation lead to 

{e} = [5]{v'} and {sCw)} -[B] {<) 
where, 


(3.33) 


[B] = 


a- 

0 

^2 

a 

0 

cN^ 

a 

0 

0 

dN^ 

dz 

0 

a 

0 

dN^ 

a 

1 ai, 

1 dN^ 

1 ail 

1 ^2 

1 ^9 

1 <^9 

42 dz 

42 a 

42 a 

42 a 

^/2 a 

42 a 


0 


0 - — 


0 

r 


r 


r 


'0 ^zz 

= [11 

0 1 ] } 

= {mf[B]{v^} 




where. 


{m} = 


(3.34) 

(3.35) 

(3.36) 


(3.37) 


Similarly, 

s,,(w) + Sgg(w) + e^{w) = { [5 ] { } 

Substituting the above equations into the eqn (3.18) leads to the following finite 
element equation in the local variables • 
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MAE NAE ( \r 

s - SK‘r{/;} + IK}’’{/=‘)+ 2Kf{/> (3.38) 

P=\ /}=! 


where, {w^} 






(3.39) 


(3.40) 



0 IK] 
[k%] [kl\ 


(3.41) 


[iy=-J (N,){mVm^’^rdrdz (3.42) 
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{w}^= (1 10 1} (3.43) 

[C] = J M{Bf{B]2 7vrdrdz (3.44) 

A‘ 

[Cl =- J {Bf {m){N^Y rdrdz (3.45) 

A* 

{/.*} = (.• } 2 ^ r dr (3 .46) 

<//) = j{C}(C}"(',‘} 2;rrdr (3.47) 

r* 

{/} = \[BY'SN\v^^^{A+v^^^\cA27crdrdz (3.48) 

{ dz dr ] 


where NAE is the number of area elements and NBE is the number of boundary 
elements. 

The variables of the area integrals (3.46 to 3.47) are transformed to the natural 
coordinates (^, 7) using the following transformation; 

+i +1 

\{..)drdz=^ ^{....)\j\d^dri (3.49) 

-! -1 

where | J| is the determinant of the elemental Jacobian matrix and is given by 



Modeling of Material Removal and Surface Roughness 


51 



dz 

dz 


drj 

dr 

dr 


drj 


(3 50) 


The choice of local coordinate system is dictated by the Gauss quadrature rule [35] 
used in the numerical evaluation of integrals over the element. 

Similarly the boundary integrals (eqns 3 46 and 3.47) are transformed to integrals 
over the interval (-1,1) by the relation 


r -7 


(3.51) 


where is the Jacobian for the boundary element and is given by 



dg dg 


2 


(3.52) 


The final finite element equation is obtained by assembling the elemental area and 
boundary matrices into global area and boundary matrices respectively. The elemental 

area matrices [k^] are calculated by evaluating the integrals by Guass quadrature 
integration technique. The assembly of the elemental area matrices into global area 
matrices is done by transferring the elements corresponding to a local degree of freedom 
in each elemental area matrix to positions of corresponding global degree of freedom in 
the global area matrix. The assembled finite element equation may be written as 

[K]{Y} = {i^}+{FJ (3.53) 


where [K] is the global coefficient matrix, Y is the global nodal primary variables 
(pressures and velocities) and F and are the global right side vectors. Here, F^ denotes 


the elastic force component. 

After substitution of boundary conditions as discussed in the next section, the final 
equation (3.53) is solved by Householder Method [36]. 


3.2.6 Boundary conditions 

In the present work, axi-symmetric abrasive media flow is considered for the 
analysis taking into consideration the schematic diagram of the AFM process shown in 

.:SNTRAL UBRAQV 

I. I. T.. KAW Wja 



^ ,1. HJLO 
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Fig. 1.4. The domain along with the boundary conditions is shown in Fig. 3.1 and 
boundary conditions are described as given below: 

(i) Entry (atAB) 

Vr = 0 ; Vz =Vp 

where Vp is the velocity of piston or velocity of flow at inlet. 

(ii) Exit ( at HG) 

Vr = 0 ; tz = 0 

(iii) Surface of media cylinder (at BC and FG) 

Vr = 0 ; t 2 = 0 

(iv) Media tooling interface (CD and EF) 

The normal component of velocity (vn) at any point on the boundary must be zero. 

Hence, 

Vr + Vz tan a = 0 ; at CD 
Vr - Vz tan a = 0 ; at EF 

where a is the angle made by the outward unit normal at tool surface with the r -axis 
(Fig. 3.1). The secondary boundary conditions are specified in terms of a friction 
condition. The Coulomb friction relation is used here. Thus 

where t^ and t„ are the tangential and normal stress components along the boundary CD 

and EF (Fig. 3.1 and 3.2). / is the coefficient of friction assumed to be constant along 
the interface. 

(iv) Surface of workpiece (DE) 

Vr = 0 ; itzi=/|tr| 

(vi) Axis of symmetry (AH) 

Vr = 0 ; tz = 0 

The geometry of flowing media is discretized using a mesh of 5x20, 9 node 
elements in the r and z directions, respectively (Fig. 3.3). The mesh emphasizes radial and 
axial refinement near the singularity. The experience with inelastic flow [30] has shown 
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that radial and axial mesh refinement is very important to resolve adequately steep 
gradients Also it has been known in viscous calculations that elements near the singularity 
with large aspect ratios produce numerical instabilities. The fine mesh offers a 
compromise between these two concerns. 


r 



Fig. 3.1 The domain and boundary conditions for axisymmetric flow of media in 

AFM process 


K 



Fig 3.2 The boundary tractions at typical boundary node K 
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Fig. 3.3 The finite element mesh and a typical element ‘e’ 


3.2.7 Solution procedure 

The initial stress rate method described by Dawson and Thompson [37] is used to 
solve the present viscoelastic flow problem. A viscous approximation for steady state 
viscoplastic analysis is extended to include the effect of the elastic strains. The solution 
procedure follows with the subscript i designating the rth iteration. 

Step 1. The approximations for stress rate are initialized to zero within each element, ie 

= 0. This in turn will require the elastic forces, for the system to be equal 

;=0 



to zero. 

Step 2. The following equation is solved for velocities and pressures. 

[K] {yL, = {F + {Fj, 


(3.54) 


For i = 0 this will correspond to viscous solution. 
Step 3. Using the velocities from step 2, 

}i+i = [.5] {v®}i+i and T,+i = fx {i}i+i 

element. 


W 



are determined for each 
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Step 4. The gradient of the stress found from step 3 is determined and with the velocity 
from step 2, the stress rate is evaluated according to eqn (3 6) 



(3 55) 


Step 5. A new elastic force matrix is computed 

{^eLi " J W Inrdrdz (3 56) 

Step 6. If there is significant change in {F^} between the /th and (/+l)th iterations, return 
to step 2 If change is less than 10'^, convergence is obtained. 



Fig. 3.4 AFM process removes high spots of the surface [10] 

3.3 Mechanism of Material Removal in AFM 

The abrasion in AFM referred here is defined as the removal of solid material from 
a surface by the unidirectional sliding action of discrete particles of another material. AFM 
process grinds away the high spots of the area, yielding a more uniform and unidirectional 
surface [10]. Fig 3.4 shows the progressive AFM action of abrasive media on the extrusion 
die [10]. The basic mechanism of abrasion in general has been the subject of many 
investigators. Khrushchov and Bavichov [38] identified following two processes taking 
place when abrasive grains make contact with the wearing surface ; (1) The formation of 
plastically impressed grooves which did not involve material removal, and (2) the 
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separation of material particles in the form of micro chips. It is found [39,40] that the chip 
cutting and rubbing depend upon the shape of indenting particles In particular, spherical 
indentors have been found to show a change over from rubbing to at least partial chip 
formation when the indentation strain (defined as the depth of indentation divided by the 
diameter of the indentor) exceeds a certain value. Pure microcutting results in a volume 
loss in the form of chips and equal to the volume of the grooves produced. 


3.4 Modeling for Material Removal and Surface Roughness 

Abrasive flow machining is considered to be a scratching action performed by 
abrasive grains in the AFM media The normal force applied to a spherical grain will 
lead it to penetrate in the surface (Fig 3.5) [41]. The grain produces a groove on the 
workpiece surface whose section corresponds to the profile of the grain. As the grain is 
translated horizontally, it removes material from the workpiece surface The amount of 
stock removal is then equal to the total volume of the grooves produced on the workpiece 
surface by each grain in the media. The idealized classical model of abrasive wear [42] 
provides a theoretical basis for the determination of material removal by AFM process. If 
the number of active grains, their shape and depth of the groove produced are known, the 
volume of stock removal can be calculated 

3.4.1 Assumptions 

The following assumptions are adopted for the analysis of material removal by 
AFM process: 

1. Most of the abrasive grains are blocky crystals that may be approximated to be 
spherical in shape [41] This is also reasonable considering the large negative rake 
angles presented by the grains. Observation of abrasive particles also shows that for 
the most part, they are generally rounded, though irregular in shape and are not 
composed of acute cutting edges . 

2. It is assumed that each grain consists of a single active cutting edge. If there are 
more than one cutting edge on one grain, there is no space to store the chip between 
the cutting edges. So that cutting edge can not be active any more. As a result, the 
grain acts as one cutting edge [43]. 



Modeling of Material Removal and Surface Roughness 


57 


3 To render the mathematics tractable it is necessary to assume that the load on each 
particle is constant and equal to the average load 

4 All active grains are assumed to be of the same size, which is equal to the average of 
given range of mesh size. 

5 Every abrasive grain is assumed to achieve the same penetration depth depending 
upon the applied force for a given workpiece material.. 



Fig, 3.5 Schematic diagram of a spherical abrasive grain removing material 

from the workpiece surface. 


3.4.2 Material removal 

The normal force acting on a spherical grain will cause it to penetrate the 

surface of the workpiece When this grain is translated horizontally, the plastically 
deformed zone beneath the surface will be inclined and give rise to upward flow thus 
forming a chip which is subsequently sheared from the surface. 

The indenting force F„g (normal force) on a spherical grain of diameter is 

given as 


4 


(3.57) 
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where, cr„ is normal stress acting on the grain. 

If h is the radius of the projected area of indentation V A, t is depth of 
indentation, and is hardness of workpiece material. Then, 

= (3.58) 

From the geometry of Fig 3.5, radius of the projected area ' of the indent 
made and depth of indentation t can be obtained as given below. 

(3 59) 




(f*=) 


(3.60) 


Substituting the value of b from eqn (3.58), into eqn. (3.60), we obtain 




k 

_ 

1 4 



(3.61) 


A' 


The cross sectional area A ' of the groove generated (shaded portion of the grain, 
as shown in Fig. 3.5) can be derived from the geometry of the figure. 

d., 

, (3.62) 


K • -1 

—Sin 

4 




Therefore, volume of the material removed ( F'g ) by an abrasive grain is equal to 
the product of area A ' and length of contact L, of grain with workpiece surface. 


dl . 2.^Hdpr) „ 

=[-sm 


(3.63) 


In practice, the total material removal is made up of a number of similar 
processes. Let n be the number of abrasive grains simultaneously acting per unit area of 
contact. Then, total no. of abrasive grains (nj indenting in the workpiece surface per 
stroke are given by 
1 . 

( 3 . 64 ) 


Vp 


where, /, is stroke length, Vp is velocity of piston, Vf is velocity of flow across 
workpiece surface, and is radius of cylindrical workpiece (Fig. 3.1) 
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From continuity equation 



where, is radius of media cylinder. 


(3.65) 


By substituting the value of — from eqn. (3.65) into eqn. (3 64), can be 
calculated as given below 

= 2 n (3.66) 

Hence, volumetric material removal in /th stroke (P^ ) in AFM is given as 

V,=2 n [-Jsin-' ^t{d^-t){-f-t)] Z, (3.67) 

This simple model (eqn 3.67) suggests that material removal per unit length of 
sliding depends on depth of indentation (given by eqn.3.61), size of abrasive grain and 
active grain density for given size of workpiece and media cylinder. 



Fig. 3.6 Simplified surface geometry 
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3.4.3 Surface roughness 

AFM removes the high spots of the surface profile of workpiece, yielding a more 
uniform surface [10], It is assumed that the surfaces of the workpiece have uniform profile 
without statistical distribution, that they have an initial surface roughness R ° , and that the 
abrasives move in the length direction of the scratches [44]. There is a constant ratio 
between the peak to valley surface roughness {Ry) and average roughness (i?a) Let R‘^ 
be the surface roughness obtained after /th stroke during AFM From Fig. 3.6, the actual 
contact length (Z, ) between the workpiece and a spherical abrasive grain is given as 


V/ 

A = (-)A 








(3.68) 


where, v and V /g are shown in Fig. 3.6 and 26^ is mean angle of surface asperity. 

Substituting the value of L, from eqn (3.68) into eqn.(3.67), the volumetric 
material removal in zth stroke is given as 


^2 d] , . 




(3.69) 


Let D =2 TV nL — L 


A . 

(f” 




- M-'Xy-Ol 


(3.70) 


Hence, V=D[l-^] (3 71) 

K 

Total volumetric material removal in N number of cycles is evaluated as 

2N 

1=1 


V = 2 TV nL-^L 


4 


sm 


2jt(d^ - t) d 2// 

-A_l [1- (3.72) 




K 


Weight of matrial removed in N number of cycles is given as 
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W,„=2np. 




W '• ‘^5 \. J. 

K.. 4 


, d 2N o< 

V'(d.-0{-f-0)I [1-§1 

^ I=l K„ 


(3.73) 


where, is the density of workpiece material. 

For a machined surface, the ratio of peak-to-valley roughness {Ry) and arithmetic 


average roughness (i?^) is closer to 5 for ground surfaces and is about 10 for honed 

surfaces [45]. This ratio is assumed to be 7 for the workpiece surface used for AFM in the 
present case. The value of this ratio obtained during experimentation is also very close to 
7. Some representative values of Ry and 7?^ are given in Table 3.1 extracted from the 

Table A.3 given in Appendix A. 


Table 3.1 Representative values of R^ and Ry of workpieces 



K 

Ry 

1 

2.08 

14.49 

2 



3 

1.85 

17.18 

4 

2.08 

15.11 

5 

2.15 

17.60 

6 

1.93 

13.98 

7 

2.31 . 

13.93 

8 

1.82 

14.01 

9 

2.12 

14.98 

10 

2.01 

15.18 


The sum of volumes of all the grooves produced during AFM process is equal to 
the volumetric material removed from the workpiece surface. Hence 

Volume of the material removed in fth stroke = actual contact length x width of 
workpiece x total height of material removed 
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v,’- *2 arr,* 7 


(3.74) 


Substituting value of V, from eqn. (3.71) into eqn. (3 74), and simplifying , we find, 
D 






14;rr^/„ 

Substituting value of D from eqn.(3.70) into eqn.(3.75), 

i i [i-- 0 1 


(3 75) 


(3.76) 


A critical surface roughness ( ) may exist at the given machining conditions, 

such that the surface roughness will no longer improve beyond that. The critical surface 
roughness is equal to the depth of indentation of spherical abrasive grains in the 
machined surface. 

(3.77) 


R. 


k 

_ 

1 4 



It should be noted that when the value of R° is comparatively small, and normal 
pressure and grain size are large enough, the value of final surface roughness may be 
greater than initial surface roughness. Such experimental observations have been 
reported by Jha [24]. Hence, selection of the process parameters should be done with due 
consideration to the initial surface roughness of workpiece. 


3.5 Results and Discussion 

The simulation procedure involves the determination of stresses on the workpiece 
surface by finite element technique as described in section 3.2 of this chapter. The normal 
stress, so obtained from the flow model, is used for the estimation of material removal and 
surface finish as mentioned in section 3.4 of the chapter. The flow chart of the procedure 
followed for determination of material removal and surface roughness in AFM is shown in 
Fig. 3.7. The viscosity of media is determined experimentally and its value at various 
abrasives concentration is reported in Appendix B. The other numerical values used in the 
simulation are given in Appendix C. 
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Figure 3 8 shows a comparison of the experimental and theoretical average 
velocity of piston and the average extrusion pressure acting on the piston during its stroke 
A good agreement exists between the results from the theoretical model and the 
experimental results It clearly indicates that the proposed flow model provides a 
reasonably good estimate of the flow characteristics involved in AFM process In Figure 
3 9, the extrusion pressures are plotted against the piston velocity for various angles of 
tooling (a) (tooling directs the media to the areas of the workpiece to be abraded Fig 
1.4). It is observed that the present FEM analysis gives a linear relationship between the 
extrusion pressure and the velocity of the piston. The piston pressure increases linearly 
with piston velocity and its value is more for higher angle of tooling. Further to achieve 
the same flow velocity of media, tooling angle should be low at low pressure Thus, for a 

given piston velocity, it requires highest pressure at 90° tooling angle which is a 
commonly used configuration. Hence, the tooling should be designed to achieve the 
desired results. Normal stress also varies linearly with piston velocity (Fig. 3.10) and its 
value is higher for higher value of angle of tooling. At lower radial pressure, the 
penetration depth expected is lower. If radial pressure is very high and the abrasive grain 
size is also high, then in some cases, it may lead to the final surface finish worse than the 
initial surface finish Hence, grain size selection and design of tooling to achieve the 
appropriate radial pressure is important so that the desired surface finish can be achieved 

Figure 3.11 shows the variation in extrusion pressure values obtained for media 
subjected to various reductions. Reduction ratio expresses the reduction in cross sectional 
area of flow. It is seen that the piston pressure increases with increase in reduction ratio 
(Fig. 3.10). With higher value of reduction ratio, the rate of increase in piston pressures is 
also higher. Figure 3.12 shows the effect of reduction ratio on normal stresses on the 
workpiece surfaces for various angles of tooling. For higher value of reduction ratio, the 
normal stress at the workpiece surface is higher. This fact can explain the experimental 
observations [12] that high abrading action occurs at the cross section of high restriction 
and less abrading action in the region of less restriction (because abrasion is directly 
proportional to normal stress). The computed stresses have been used to evaluate the 
machining forces acting in abrasive flow machining. The machining force so evaluated 
from the analysis has been employed to carry out theoretical analysis of material removal 
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and surface roughness in AFM Thus, analysis of media flow is an important step in the 
direction of modeling of AFM process 



Input: 

Specifications of abrasive 
media, workpiece, tooling 
and AFM conditions 




Fig. 3.7 Flow chart for determination of material removal and surface roughness 

by theoretical model 
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Fig. 3.8 Variation in extrusion pressure with average piston velocity 
(a=30", C=66%, /?, = 90.50%) 



Piston velocity, cm/s 

Fig. 3.9 Variation in extrusion pressure with piston velocity for various 

tooling angles (a) 
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Fig. 3.10 Variation in normal stress on workpiece surface with piston velocity for 

various tooling angles (a) 



Fig. 3.11 Variation in extrusion pressure with reduction ratio for various 

tooling angles (a) 
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Reduction ratio, % 


Fig. 3.12 Variation of normal stress with reduction ratio for various tooling angles (a) 

3.5.1 Material removal 

The values of number of active abrasive grains per unit area of contact under static 
conditions at various concentrations and mesh sizes have been reported in chapter 4. The 
same values have been used in the present computations of material removal and surface 
roughness. An estimate of the depth of indentation of abrasive in the workpiece material 
reveals that its value decreases with increase in the hardness of workpiece material (Fig. 
3.13) and mesh size of grains (Fig. 3.14). Figure 3.15 compares the response surface 
analysis (RSA) and theoretical results between the extrusion pressure and material 
removal. It is found that material removal increases with increase in extrusion pressure. 
The simulation results for material removal agree well with the RSA results at low 
extrusion pressures. But at higher extrusion pressures, number of abrasive grains taking 
part in actual cutting may be less due to rolling action of grains at higher velocity. 
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Hence, RSA material removal is less than theoretical (Fig. 3.15). Due to higher pressure, 
penetration depth increases. Now, if the viscosity of media is not high enough to keep 
the grain intact with the media, due to larger tangential force it may rotate and reenter 
inside the media resulting in lower depth of indentation. Another possibility is that 
increased force may not result in increased depth of indentation. Because of low 
viscosity of media, it may not be able to withstand the acting high force, and ultimately 
the grain may penetrate inside the media resulting in the reduced depth of indentation in 
the workpiece material. Hence, lower material removal is . achieved compared to the 
theoretical value. This phenomenon of rolling and/or reentry inside the media has not 
been accounted in the theoretical analysis. Hence, the increasing difference between 
theoretical and RSA results at higher extrusion pressure is seen in Fig. 3.15 Material 
removal also increases with increase in reduction ratio (Fig. 3.16). It is due to increase in 
depth of indentation at higher reduction ratio. At higher reduction ratio, the normal stress 
acting on the workpiece surface is higher (Fig. 3.12), leading to higher material removal. 

Figure 3.17 shows the effect of percentage concentration of abrasives by weight on 
material removal. As the percentage concentration of abrasives in media increases, the 
material removal increases because higher concentration results in higher viscosity of 
media (Fig. B.3 in Appendix B). Hence, abrasive grains (or media) can sustain higher 
normal and tangential forces. As a result, higher depth of indentation (or higher material 
removal) can be achieved. With higher percentage concentration of abrasives, more 
number of abrasive grains come in contact with the surface to be machined (Fig. 4.7 in 
chapter 4). All these factors lead to an increase in material removal. However, the 
difference between RSA results and theoretical ones increases with increase in the 
percentage concentration of abrasives. It can be explained on the basis of comparatively 
low viscosity of the carrier and non-linear relationship between percentage concentration 
of abrasives and number of active grains (as discussed in chapter 4). Material removal 
also increases with number of cycles (Fig. 3.18). However, in this case also, the 
deviation between the two curves increases with the increase in number of cycles. 
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Hardness of workpiece material, N/mm^ 


Fig. 3.13 Effect of hardness of workpiece material on depth of indentation 

(/7,= 50bar,M=55) 



Fig. 3.14 Effect of average mesh size of grains on depth of indentation 
(/?^= 50 bar, = 2177.80 N/mm^) 
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Extrusion pressure, bar 


Fig. 3.15 Variation of material removal with extrusion pressure 

0.805, 15, C = 66 %) 



0.84 0.8725 0.905 0.9375 0.97 

Reduction ratio 


Fig. 3.16 Variation of material removal with reduction ratio 
50 bar, 15, C = 66 %) 
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3.17 Variation of material removal with percentage concentration of abrasives 

( /?, = 50 bar, = 0.905 N = 15) 



Fig. 3.18 Variation of material removal with number of cycles 
( /jg = 50 bar, = 0.905, C = 66 %) 
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3.5.2 Surface roughness 

Figure 3 19 exhibits the role of extrusion pressure on surface roughness value 
obtained after the specified number of machining cycles. It is observed that change in 

surface roughness value ( AR^ = R° -R'^) increases with increase in extrusion pressure In 

response surface analysis (RSA) results, increase in change in surface roughness at higher 
pressure is less than theoretical value. It may be due to rolling action of grains at high 
pressure Change in surface roughness increases with increase in percentage concentration 
of abrasives (Fig 3.20) With higher percentage concentration of abrasives, the number of 
abrasives taking part in machining will be more, and hence surface roughness obtained for 
a specified number of cycles will be better. Change in surface roughness also increases 
with reduction ratio (Fig 3 21) and number of cycles (Fig 3.22), due to increase in length 
of slug of media passed across the workpiece surface. It has been found theoretically that 
for workpiece of low hardness material, the surface roughness value decreases upto a 
certain value of extrusion pressure for specified number of cycles (Fig. 3.23). Beyond this 
value of extrusion pressure, the surface value can not be further improved. If the pressure 
is large enough, then it may lead the surface roughness worse than the initial one due to 
increase in depth of indentation 



Fig. 3.19 Variation of change in surface roughness with extrusion pressure 

(R, = 0.805, N=15,C =66 %) 
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Fig, 3,20 Variation of change in surface roughness with percentage concentration of 
abrasives (Pe= 50 bar, 7?^= 0.905 N= 15) 
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Fig 3.21 Variation of change in surface roughness with reduction ratio 
(p,^50 bar, C=66%) 
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Fig. 3.22 Variation of change in surface roughness with number of cycles 
[ /?, = 50 bar, = 0.905, C =66 %J 



Fig. 3.23 Variation of surface roughness value with extrusion pressure 

= 500N/mm\ N= 10, 7?, = 0.905, C =66 Vo, M = 50-60) 
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3.6 Conclusions 

In the present chapter, a model for the flow of AFM media through cylindrical 
workpiece has been developed and solved by finite element method It predicts the normal 
stresses at the workpiece surface A model for estimation of the material removal and 
change in surface roughness has also been presented The analysis includes the operating 
conditions of the process, viz abrasive size, number of contacting abrasive particles, 
extrusion pressure, workpiece material hardness, and size of workpiece to be machined 
Good agreement is found between the RSA results (which are derived from experimetal 
data) and theoretical predictions. The following conclusions have been derived fi'om the 
present analysis 

1 . A good agreement exists between the theoretical and experimental results for extrusion 
pressure and average piston velocity It clearly indicates that the proposed flow model 
provides a reasonably good estimate of the flow characteristics involved in AFM 
process 

2 The piston velocity increases linearly with extrusion pressure and its value is higher 
for larger angle of tooling. 

3 At higher value of reduction ratio, the rate of increase in extrusion pressure is higher 
The normal stresses on the workpiece surface increase with increase in reduction ratio 

4 Material removal increases with extrusion pressure and percentage concentration of 
abrasives in media. 

5. Change in surface roughness value increases with increase in extrusion pressure, 
percentage concentration of abrasives and reduction ratio for a specified number of 
cycles 

6. Beyond a certain value of extrusion pressure, surface roughness deteriorates due to 
increase in depth of indentation. The minimum surface roughness value is equal to the 
critical surface roughness. 



Chapter 4 

Simulation of Surface Generation in AFM 


4.1 Introduction 

The theoretical model discussed in the previous chapter has been derived, assuming 
constant depth of indentation of every grain and uniform surface profile without statistical 
distribution. Stochastic simulation methodology has been discussed in this chapter to 
generate surface profile which is more close to actual AFM results. The process of finishing 
in AFM is random removal of material by individual abrasive grains in the form of 
microscopic chips. The shape of removed material is determined by the kinematics and the 
shape of abrasive grains. It is easier to calculate the shape of cut groove from the relative 
motion between the workpiece material and single grain under consideration. However, it is 
difficult to estimate the roughness of the machined surface. This difficulty is due to the 
complicated and random distribution of grains in the media. Abrasive flow machining is 
viewed as a stochastic process whose events are generated by the interactions between the 
abrasive grains and workpiece. By simulating these individual events, it is possible to 
predict the characteristics of surface likely to be generated in AFM. Simulation techniques 
can help in understanding the mechanics of AFM process by predicting the effects of 
various parameters on the responses and comparing them with RSA results. 

Not much information is available in the literature which deals Avith the simulation 
of surface generation and material removal in AFM process. In case of grinding, there have 
been many attempts in the past to model the surface generated using statistical approach 
[46]. Yoshikawa and Sata [47] simulated the grinding process by the Monte Carlo method. 
Random numbers were used to give the coordinates of the position of a grain on the 
grinding wheel. By this method, the distribution of the grains spacing and of the 
undeformed chip geometry, together with the generated surface profile are obtainable from 
the structure of the grinding wheel which has random character. The assumption of the 
grains existing independent of each other on the wheel surface is the major drawback of this 
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model. A grinding model is postulated by Law and Wu [48], in terms of grain distribution 
on the abrasive wheel and kinematic grinding conditions. Transverse workpiece profiles 
have been conjectured as sample realizations of the stochastic grinding process. Hamed et 
al. [49] applied an approach to generate random surfaces by computer simulation which 
involves the use of a series of 'unit events' to produce a surface profile. The work 
confirmed that the random engineering surfaces can be described statistically by a profile 
having ordinate heights with a Gaussian distribution together with a form of exponential 
autocorrelation function Bhateja [50] evolved a hypothesis that the ground surface texture 
is complementary to the envelope of a large number of active transverse wheel profiles that 
traverse the grinding zone. Based on this approach, the relationship between the roughness 
of the wheel and the workpiece could be established. Classical two body abrasion theories 
have been applied in a statistical simulation model by Staffan et al. [5 1]. A unique feature of 
the model is that the abrasive simulation is performed on a surface of realistic topography, 
developed during the simulation process. Koshy et al [46] presented a stochastic simulation 
of the diamond grinding wheel which entails a random sampling technique. The indices that 
characterize the topography of a diamond wheel viz. the static grain count, the protrusion 
height distribution, the distribution of inter grains spacing, and the projected area due to 
exposed abrasives are estimated by this model. 

The approaches discussed above, to model surface generation can be applied to 
AFM process with certain modifications depending upon the AFM conditions. As the 
distribution of grain cutting edges in the media are random and grains which contribute to 
generate a measurable profile on the workpiece are numerous, it is very difficult to describe 
the necessary configuration of media surface in a deterministic way. On the other hand, it is 
possible to obtain a complete profile of workpiece surface generated by motion of the 
grains, by simulation. The surface generated in AFM is the cumulative outcome of 
numerous abrasive-workpiece interaction events. By simulating these individual events, it is 
hence possible to predict the characteristics of machined surface. Implementation of the 
same, requires information on the topography of the AFM media and the kinematics of the 
AFM process. 

The present chapter describes the concepts of a proposed versatile simulation 
model, which generates and statistically evaluates the interaction between spherical abrasive 
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grains and a workpiece surface The simulation enables prediction of the roughness of the 
machined surface and material removal with reference to the abrasive grit size and 
percentage concentration, in addition to obtaining the statistical characteristics of machined 
surface. Active grain density has also been estimated from the present stochastic simulation. 


4.2 Simulation of Surface Generation 

The method applied here to explore the nature of the generated profile is simulation 
by employing mechanics of AFM process. The model used for material removal in AFM 
process is restricted to pure microcutting without material displacement. All considerations 
are geometrical in nature. During AFM process, media is passed across the workpiece 
surface. AFM is considered to be a scratching action performed by abrasive grains in the 
media. The force applied to the media will cause an active grain to indent in the workpiece 
surface. The grain produces a groove on the workpiece surface and the groove section 
corresponds to the profile of the grain. As the grain is translated horizontally, it removes 
material from the workpiece surface. 

The simulation of surface generation in AFM is carried out based upon the 
following assumptions: 

1 . The path traced by an individual grain is straight line. 

2. The material is removed in the form of chips without side pile up, whenever grain 
workpiece interference occurs. 

3. Abrasive grain is approximated by a sphere. This is also reasonable considering the 
large negative rake angles presented by the grains [41]. Observation of abrasive 
particles also shows that they are generally rounded in shape, and are not composed of 
acute cutting edges [52]. 

4. During the manufacturing of particles, a large number of chance factors are at work. 
This leads to the emergence of a normal distribution of particle size, which has been 
observed to hold for narrow range [46]. Hence, the distribution of abrasives’ radii is 
assumed to be normal, and symmetric about the mean grain radius. 
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Since each nominal grit size includes a range of abrasive particle size (Fig. 4.1), two 
grit numbers and A/j, are used to denote the abrasive mesh size. refers to the 

sieve that would let the abrasive grains pass through and refers to the sieve that would 

detain most of them The grit radius is related to the sieve number [3] as 

Tg (mm) = 7 62 

As the distribution of abrasives radii is assumed to be normal and symmetric about 
its mean, the mean grain radius is given as follows. 

V - ^so - 3 31 (4 2) 

2 MM, 



Fig. 4.1 Schematic representation of distribution of grains in abrasive media 


Since, the range of normally distributed random variables, in general, is six times the 
standard deviation, the standard deviation of the grain radius is evaluated as 


1.27 


(4.3) 


MM, 

For AFM media having C % concentration of abrasives by weight, the total volume 
of abrasive grains (F^ ) in the volume of media is given by 






Fig. 4.2 (a) Schematic representation of simulated AFM media 

To simplify the analysis, total periphery of circular cross section of the workpiece is 
converted to equivalent periphery of square cross section (Fig. 4.2 (a)); keeping total length 
of media passed through the workpiece as constant. 


Side of square cross section, w 


2 


(4.7) 
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Fig. 4.2 (b) Cross sectional representation of workpiece surface and media 


Let a square cross section of side w and length 4 (Fig.4 2) of the AFM media is 
considered for simulation. Then, volume of abrasive grain in the volume of abrasive media 
(^m~ /„) passed in one stroke across workpiece surface is related as 


K = 


loop. 


(4.8 a) 




C LPc 


Cp,+im-c)p, 


(4.8 b) 


Only one side out of the four sides of a square cross section has been considered 
for simulation It is assumed that there is same surface pattern and amount of material 
removal obtained for the remaining three sides of workpiece. This saved the computational 
time with a small difference in the results obtained. During computation, it is found that 
total number of grains in the media passed in one stroke across the workpiece surfaces are 
of the order of for the conditions taken in the AFM simulation. It creates 
computational and memory problems. Hence, a thickness of 5 mm of abrasive media has 
been considered to be flowing over workpiece surface, keeping the length and width of 
flowing media constant. This reduction in simulated thickness of media does not affect the 
number of grams interacting with surface profile and distribution of grains in the media near 
the workpiece surface. On the other hand, this modification reduced the computational time 
and memory requirements drastically 
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The first step in the simulation is to generate 'N^' normally distributed random 
numbers with mean radius and standard deviation which correspond to the grain 
radius such that 

4 Na 

3 ^ ,5 '■g: K (4 9) 

where is the radius of the /th grain and is the total number of grains in the simulated 
volume 


For a spherical grain of diameter , the maximum depth of indentation t of a grain 


in workpiece material has been derived in chapter 3 and can be given by eqn (4. 10). 



(4.10) 


where a„ is the normal stress acting on the grain and is the hardness of workpiece 
material. 

The distribution of grains in the media is represented by giving their coordinates (x, 
y, z), the components of which give the position of a grain in axial, transverse and 
perpendicular to surface directions. Grains are randomly positioned in the AFM media so 
that X, y, and z generally have a uniform distribution. Thus, the position of the grains with 
respect to the surface is decided by a random number generator. A subroutine from NAG 
software is used for the generation of random numbers. 

Now, a set of these uniformly distributed random numbers N^ (xj , yi , Zi ), is 


generated to represent the x,y,z coordinates of the center of each of N^ grains. 

Only grains which do not overlap are considered. For this purpose, the grains are 
assumed to occupy a region of influence equivalent to nominal grain diameter. Generation 
of the coordinates of the fth grain is therefore subject to the condition, that distance 
between two adjacent grains is equal to or more than the sum of their radii. 

[(>‘,-x,)"+(y,-y,)"+(z,-z,)"]“2(rj,+r^) ; y=lto/-l (4.11) 

The profile of the workpiece surface can be displayed in the form of a two 
dimensional cross-section perpendicular to the grooving direction. In the manufacturing 
processes involving the random removal of material, the roughness is effectively produced 
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by many independent particles’ hits Hence, the amplitude distribution of the roughness 
waveform could be justifiably expected to be Gaussian no matter how much of the 
individual hits have affected the surface profile [48], The original profile of the workpiece 
is taken to be Gaussian, represented in digital form and coordinates (x',z') of the profile are 
stored in the computer at 0 1 /an spatial resolution. 

The centers of the abrasive grains (x,, 2 ,) are defined with respect to the coordinate 

system fixed to the workpiece. Active grains are only those grains, which actually take part 
in the material removal In Fig 4.3, grain G1 is interacting with workpiece profile, hence it 
is an active grain Whereas, grain G2 is passing over the workpiece surface profile without 
removing any material It is considered as an inactive grain. During the simulation process, 
if the active grain interacts with the workpiece profile at the appropriate location, then a hit 
is recorded. The grain indenting the workpiece surface would have the z coordinates of its 
surface below the corresponding z coordinates of surface profile of the workpiece. These 
grains can be identified by segregating grains which satisfy the following condition. 

z' - (z, -r^,) > 0.0 (4.12) 

The machining process is simulated by comparing the coordinates of a grain with 
the corresponding surface points on the workpiece. The coordinates of the cutting surface 
of a grain are generated by solving the equation of a circle for discrete values of x. An 
active grain such as G1 of radius r, (Fig.4.3) would modify the work profile only within the 

interval (x, - r, ) to (x, +rj . It is assumed that the cross-sectional shape of the groove 

formed by the sliding abrasives on the workpiece surface is truncated circular profile. In 
other words, material removal is by ideal ductile microcutting with no pile-up of material 
on the sides of the groove. 

Material removal in AFM process is simulated by successive abrasive grain 
grooving in the workpiece. The height of the profile at the location is then modified to 
correspond to this interaction. The modification of the profile of the work surface as each 
active grain slides past the work surface is modeled by transforming the coordinates of the 
profile of individual active abrasive grain to the workpiece, based on geometrical 
considerations alone. When a grain passes across the material, the computer calculates the 
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Fig. 4.3 Schematic diagram of interaction of an abrasive grain with workpiece 

surface 

successive new heights of surface profile from the reference surface, thus generating new 
surface profiles. After each grooving event the groove width and volume of material 
removed are calculated and stored. After completion of simulation of abrasion by each 
abrasive grain, the resulting surface profile is calculated. Thus, a simulated sample profile of 
the machined surface, and volume of material removed after each stroke are obtained. The 
whole process of simulation from generation of random numbers to modification of 
workpiece surface profile is continued until the required number of cycles of AFM process 
is completed. The flow chart of AFM simulation is shown in Fig. 4 4 

For studying the characteristics of the machined surface, a projected profile length of 
12 mm was considered at a resolution of 1 fim. For sampling length, the mean line was 
computed such as the area of the solids above the line is equal to the area of the voids 
below it The roughness parameter is the average height from a mean line of all 
ordinates of the surface, and is calculated as follows; 

1 ^ 

L- 0 

where z(x) is the height of the profile from the mean line and L is the sampling length 
used to compute the average roughness of the profile. Since the initial surface roughness of 
all the workpieces was not exactly same, the change in surface roughness of the machined 
surface has been plotted for comparison purposes. 
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Fig. 4.4 Flow chart of AFM simulation 
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4.3 Active Grain Density in AFM Media 

The spatial distribution of abrasive grains over the media surface in contact with the 
workpiece and their morphology comprise the AFM media topography The AFM media 
topography is characterized microscopically by nominal grain size, the static grain density 
and a three dimensional distribution of grain positions. 

The cutting edges in media are geometrically undefined in location and shape As a 
result, estimations are difficult to make with any degree of accuracy It has, therefore, been 
necessary to develop techniques for examining details of the media topography in order to 
ascertain quantitative measurements of cutting edges. 

A method has been proposed by Rajurkar [17] to estimate the number of abrasive 
grains by using the frequency associated with the primary root from the data dependent 
system (DDS) models of the AFM generated surface profiles. The range of values of 
number of dynamic active grains (for 70 mesh size with 66 % concentration of abrasives in 
medium viscosity media) over a unit area of extrusion passage (in^) was found to be in the 
range of 98 to 3 12 Values of active grains in media at other concentrations and mesh sizes 
have not been reported. No systematic effort has been made to estimate the topography of 
media and simulate the generated profile theoretically. 

But characterization of the topography of conventional grinding wheels has 
introduced considerable interest among researchers in the past Reference [46] presented a 
survey of models available for the topography of conventional grinding wheel. Nakayama 
and Shaw [53] proposed the distribution of effective grain density versus the radial distance 
from the outermost grains by a dynamic method. In these studies, the roughness equations 
derived or grain distribution employed contain parameters which are difficult to determine 
by experiments. 

The objective of this section is to present a method, utilizing microscope techniques 
for quantitatively characterizing the topography of AFM media. The stochastic simulation 
method described in the previous section is also used to evaluate the active grain density. 

Scanning electron microscope (SEM) has been tried for the analysis of topography 
of media. It is observed that the shape of a sample surface of polyborosiloxane media keeps 
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changing under vacuum which is applied in SEM It is due to the flexible and semisolid 
nature of polyborosiloxane media Hence, it is not possible to analyze correctly the 
topography of media by SEM. As a solution to this, optical microscope is used for 
determining the active grain density of media 

4.3.1 Active grain density measurement by optical microscope 

The undefined cutting area and the extremely small dimensions of the abrasive 
grains under observation present barriers in the assessment of the cutting surface of media. 
It is therefore suggested to use a microscope for the analysis of the cutting surface of the 
media 


During AFM only a small number of grains come in contact with the work material. 
Whether a grain will be active or not, depends on the kinematics and geometric properties 
of the process. Hence, active grains are those which are supposed to take part in material 
removal process or those particles contacting the workpiece surface during machining. The 
static grain count is the number of all protruding grains on the media surface per unit area. 

Microscopic method involves the use of an optical microscope or scanning electron 
microscope to observe and measure topographic features of the media surface. When 
viewing normal to the media surface, flattened areas on the grain tips reflect light and 
appear shiny against a dark background [3]. Such areas are often referred to as wear flats 
since during movement of media across the surface to be machined, wear of grains by 
rubbing takes place. 

The media with different grit sizes and concentration of abrasives have been 
employed to determine the number of cutting edges on the media surface. The media 
leaving the workpiece surface after machining is taken out of the AFM set-up without 
disturbing the cross sectional shape and position of grains, and then the value of grain 
density is measured by means of an optical microscope sighting on the distinct grains. Since 
the microscopic measurements are made on media surface after machining, the 
concentration of the cutting points (or grain densities) and their spacings represent dynamic 
values for the AFM conditions. 
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For the purpose of measurement and analysis, a number of viewing areas is 
randomly selected. The number of distinct grains counted per unit area of viewing on the 
media surface is interpreted as “active grain density”. The average number of cutting points 
per unit length can be estimated as the number of grains per unit length intersecting a 
straight line superimposed in the viewing area 

The values of active grain density for various concentration and mesh size, for two 
different carriers (putty and polyborosiloxane), are given in Table A 5 (Appendix A) The 
abrasive used with both carriers is silicon carbide 

4.3.2 Results and comparison with stochastic simulation 

Fig 4.5 represents the variation in active grain density as a function of abrasive 

concentration for different grit sizes wdth putty as carrier. The active grain density increases 

with an increase in abrasive concentration, and for a particular concentration, it is higher for 

a larger mesh number (or smaller grain diameter) (Fig. 4.6). Fig. 4 7 shows the variation of 

active grain density with concentration when polyborosiloxane is used as a carrier. Only 

mesh size of 50-60 is used with polyborosiloxane due to limited availability of abrasive 

media with polyborosiloxane which is procured from Exrtrude Hone, USA. The active 
grain density at a particular percentage concentration of abrasives by weight is higher for 

putty as carrier than for polyborosiloxane (Table A. 5, Appendix A). As density of putty 

(2240.0 kg/m^) is higher than polyborosiloxane (985.0 kg/m^), for the same weight of 

media, the volume of polyborosiloxane carrier is higher than putty. As a result, the number 

of abrasive grains per unit area will be lower in polyborosiloxane abrasive mixture than 

putty abrasive mixture at a particular percentage concentration of abrasives by weight. As 

seen from the Figs. 4. 5-4. 7, the experimental values are generally lower than the calculated 

values of active grain density. One reason for this may be experimental difficulties in 

detecting all contacts and an unintentional tendency to overlook small contact areas when 

scanning large abrasives. However, the simulated results and the results obtained from the 

microscopic measurement method are close to each other, and hence validate the simulated 

results. 
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In actual practice, the value of active grains should remain relatively stable for low 
number of strokes. There may be slight increase in the active grain density, as grains 
fracture produces more smaller grains Once the number of strokes reaches a critical value 
and the media has absorbed too much of chips, the value of active grain density may take a 
drastic downturn. Hence, media must be replaced after a certain number of cycles to 
maintain the performance of the process Further work is required to determine the grain 
wear and critical number of cycles for media replacement. 



Fig. 4.5 Variation of active grain density with percentage concentration and mesh 

size for putty as carrier 
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Fig. 4.6 Variation of active grain density with average mesh size for putty as carrier 

(C=40%) 



Fig.4.7 Variation of active grain density with percentage concentration of abrasives 

for polyborosiloxane as carrier 
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4.4 Results and Discussion 

The conditions, under which the present simulation has been earned out, are given 
in Appendix C. The assumed initial Guassian surface profile of the workpiece is shown in 
Fig 4.8(a) The surface profile obtained after simulation fot the specified machining 
conditions is illustrated in Fig 4 8(b). The sampling length is 1 mm with coordinates at 1 
fm resolution This profile resembles with the actual surface profile obtained from AFM 
experiments. Considerations of real shape of grains, eftect of weai of giain, eftect of 
elasticity, pile up of material and actual initial surface profile (.if the work can piovide more 
realistic results in the simulation in the succeeding sections, the lesults foi material 
removal and surface roughness, predicted by simulatitm, are discussed and compared with 
the response surface analysis (RSA) results and theoretical icsults (discussed in chajiter 3) 

4.4.1 Material removal 

Figure 4.9 shows the eftect of percentage camcentiation of alirasives by weight on 
material removal As the percentage concentration of abrasives in media inci eases material 
removal also increases, because higher concentration of abrasives lesults in higher number 

of active grains and hence higher material removal 'I'he diftercncc lictween the RSA results 
and results due to theoretical model and simulation may tie due to inaccuracy in 

measurement of active grain density specially at higher value of concentration Also 

generally, the measured value and stochastically determined value of active grain density is 

higher than the effective active grain density during the AFM process I Icnce, tfie predicted 

values of material removal at higher concentration are generallv highci than the RSA data, 

Figure 4.10 shows the effect of mesh size on material removal for mild steel 
workpieces. From the figure, it is evident that as the mesh size increases, the material 
removal increases. The reason for increase of material removal is that with increase in mesh 
size, area of the indented grooves and depth of indentation decrease but total number of 
interferences between grains and workpiece increases due to increase in active grain density 
(Table A.5 Appendix A and Fig 4.6). The overall eftect is an increase in material removal 
with increase in abrasive mesh size. 
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It is found that material removal increases with increase in extrusion pressure (Fig. 

4 11) The simulation results for material removal agree well with the RSA data. But at 
higher pressure, the number of abrasive grains taking part in actual cutting may be less due 
to rolling of grains at higher velocity. Hence, actual material removal will be less than 
theoretical one which can also be seen from Fig 4. 1 1 also. This deviation would increase 
with increase in extrusion pressure Material removal also increases with increase in 
reduction ratio (Fig 4 1 2), due to increase in the slug length of media passed per unit time. 
At higher reduction ratio, the stresses acting on the workpiece surface are higher, leading 
to higher material removal. At lower value of reduction ratio, the predicted material 
removal can be observed to be in good agreement with the RSA results. At higher end of 
reduction ratio, a difference between the results exists. At higher reduction ratio, more 
grains come in contact with the workpiece surface in a cycle, due to increase in slug length 
Secondly, normal stresses at workpiece surface are higher at higher reduction ratio for a 
specified extrusion pressure. As a result, material removal is higher due to higher depth of 
indentation at higher normal stress. The same is also true with higher extrusion pressure. 
Material removal value increases with the number of cycles (Fig.4 13) and also the 
difference between the simulated and RSA results increases with increase in number of 

cycles. The difference in RSA values and theoretical results is more at higher number of 
cycles. Material removal in a specified number of cycles decreases with increase in hardness 

of workpiece material due to decrease in depth of indentation (Fig. 4.14). 

4.4.2 Surface roughness 

The effect of percentage concentration of abrasives on the change in surface 
roughness indicates increase of change in surface roughness (ARJ with increase m 
percentage concentration of abrasives (Fig. 4.15). With higher concentration of abrasives, 
the number of abrasives taking part in machining will be more, and hence surface roughness 
obtained for a specified number of cycles will be better. The values due to theoretical model 
are higher than that predicted by simulation for all concentrations. The predicted AR^ can 
be observed to be in better agreement at low percentage concentrations. At higher 
percentage of abrasives concentration, the AR^ is higher. This may be due to higher value 
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Fig. 4.9 Variation of material removal with percentage concentration of abrasives 
(/?,= 50 bar, M =50-60, R^= 0.905, V=15) 



Fig. 4.10 Simulated variation of material removal with mesh size of grains 
(/?,= 50 bar, = 0.905, V=15, C= 66%) 
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Fig. 4.11 Variation of material removal with extrusion pressure 
(M =50-60, 0.905, A^=15, C=66%) 



Fig. 4.12 Variation of material removal with reduction ratio 
(/7g= 50 bar, M =50-60, A^=15, C = 66%) 
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Fig. 4.13 Variation of material removal with number of cycles 
( /?, = 50 bar, M =50-60, i?, = 0.905, C = 66%) 



Fig. 4.14 Simulated variation of material removal with hardness of workpiece 

material 

(;?^=50 bar, M =50-60, i?,= 0.905, iV=15, C=66%) 
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of active grain density predicted by stochastic simulation and those measured by 
microscopic method, at higher percentage concentration. 

AR^ value is higher for increased mesh size (Fig 4 16) If the grain size is large 
enough, then it may lead to the surface roughness worse than the initial one due to large 
depth of indentation. This can be explained with the help of simulation model also 

Figure 4 1 7 exhibits the role of extrusion pressure on change in surface roughness 
value obtained after the specified number of machining cycles. It is observed that AR^ 

value increases with increase in extrusion pressure Figure 4.17 also presents a comparison 
of the change in surface roughness as a function of extrusion pressure, computed by 
theoretical model, simulation and due to experiments. The simulated results are in 
agreement with those of the RSA results. However, difference is more in case of the 
theoretical model This may be due to the assumption taken in modeling that only surface 
peaks are removed during AFM, but at higher pressure, in some sections of peaks, there 
may not be uniform material removal due to higher depth of indentation. As a result there is 
less change in surface roughness than those obtained by theoretical model. AR^ value also 

increases with number of cycles (Fig 4.18) and reduction ratio (Fig. 4.19). Comparison of 
results indicates that results due to theoretical model, simulation and RSA are in good 
agreement as far as lower numbers of cycles are concerned At higher value of number of 
cycles, difference between the RSA results and results obtained from the theoretical model 
and simulation increases. In practice, there may be more rolling of grains than sliding with 
reduction in surface roughness value at higher number of cycles. As a result, there will be 
less change in surface roughness at higher number of AFM cycles. Change in surface 
roughness decreases with increase in the hardness of workpiece material due to decrease in 
depth of indentation (Fig. 4.20). 
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Fig, 4.15 Variation of change in surface roughness with percentage concentration of 
abrasives (/?g= 50 bar, M =50-60, 7?^= 0.905, N =15) 



Fig. 4.16 Simulated variation of change in surface roughness with mesh size of 

grains 

(Pe==50 bar, = 0.905, N=15, C= 66%) 
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Fig. 4.17 Variation of change in surface roughness with extrusion pressure 
(M =50-60, i?,= 0.905, N=15, C=66%) 
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Fig. 4.18 Variation of change in surface roughness with number of cycles 
50 bar, M =50-60, 0.905, C= 66%) 
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Fig. 4.19 Variation of change in surface roughness with reduction ratio 
(/?^= 50 bar, M =50-60, iV= 15, C= 66%) 
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Fig. 4.20 Simulated variation of change in surface roughness with hardness of 
workpiece material {pg=S0 bar, M =50-60, i?g= 0.905, C= 66%, N= 15) 
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4.5 Conclusions 

The generation of surface profile and material removal in AFM are simulated 
considering random distribution of abrasive grains in media and their interaction with 
workpiece profile during the process. The results of the proposed simulation and RSA are 
in good agreement, which justifies the use of this approach for generating the surface in 
AFM The simulation enables the prediction of surface roughness and material removal 
with reference to percentage concentration and mesh size of abrasives, extrusion pressure, 
number of cycles and reduction ratio The active number of abrasive grains per unit volume 
of media flow increases with percentage concentration of abrasives and reduction ratio, 
hence higher material removal and change in surface roughness is achieved Extrusion 
pressure, grain mesh size and hardness of workpiece affect the depth of indentation of a 
grain in workpiece material The significance of these parameters on final surface profile is 
discussed with reference to AFM process Change in surface roughness value and material 
removal increases with extrusion pressure, number of cycles, reduction ratio and 
percentage concentration of abrasives 



Chapter 5 

Thermal Aspects of AFM Process 


5.1 Introduction 

The AFM process modeling is complex due to little understood behaviour of the 
media and the complicated and random nature of material removal Both theoretical and 
empirical studies of abrasive flow machining are greatly hampered due to the inherent 
random nature and multiplicity of variables. To date it appears that there are no theoretical 
modeling for determination of tangential forces, specific energy and thermal aspects of 
AFM process. There are no systematic experimental measurement of forces and 
temperature distribution during AFM and its dependence on the relevant process 
parameters Thus, studies aimed at predicting the AFM performance are greatly hampered. 

Inspite of the process being a critical one for the industries, the mechanics of the 
AFM process has not been studied in any great detail. In this study, an expression has been 
derived for the determination of tangential forces and specific energy. The tangential force 
acting during AFM is explained in terms of an equation which takes into account the 
forces arising from material removal and from rubbing between the workpiece and 
abrasives. The understanding of the forces provides a basis for the investigation of the 
AFM process. Preliminary results are also presented for the effect of process parameters 
on force ratio and specific energy. These results will allow the mechanics of AFM to be 
compared to the mechanics of other abrasive processes such as grinding. 

An investigation on the theoretical and experimental determination of thermal 
properties of media has been carried out by Fletcher [22]. Davies and Fletcher [20] found 
experimentally that viscosity of media is significantly affected by temperature. Viscosity 
of high viscosity media reduces drastically even with a small increase in temperature (2-10 
°C). The results [20] have shown that progression from low to high viscosity base media 
produces a reduction in the temperature rise. Hence, it is of great practical and theoretical 
interest to study the temperature rise during the AFM process. In the case of AFM, the 
complex nature of the process precludes a precise thermal analysis, and consideration of 
simplified models is required. The purpose of this chapter is to present a simple theoretical 
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thermal analysis and to provide a means for predicting approximate temperatures during 
the AFM process. 


5.2 Forces and Specific Energy in AFM 

In this section the mechanics of AFM process is analyzed from the consideration of 
the forces acting on a single grain in machining. A simple analysis based on abrasion 
theory for calculating the cutting forces due to a single abrasive grain can be expanded to 
predict the forces arising in AFM The machining action during AFM compares to a 
grinding operation as media uniformly removes material from the workpiece surface. 
When media is forced through a restricted passage, the media behaves as a deformable 
grinding stone [16]. Shaw [41] proposed a new theory of microchip formation for fine 
finish grinding assuming spherical abrasive grain. Brian Rowe and Xun Chen [54] also 
applied this theory for analysis and simulation of grinding process. In abrasive flow 
machining process, the normal force {F„^) applied to such a spherical grain will cause it to 

penetrate the surface (Fig. 5 1). As the grain is translated horizontally, the plastically 
deformed zone beneath the surface will be inclined and will give rise to upward flow, thus 
forming a microchip which is subsequently sheared from the surface The following 
assumptions are made for determining the specific energy in AFM; 

1 . The path traced by an individual grain is straight line 

2 The material is removed through plastic deformation in the form of chips without the 
side pile up, whenever grain workpiece interference occurs 
3. Abrasive grain is approximated by a sphere [41,54], This is also reasonable, 
considering the large negative rake angles presented by the grains. 

Fig. 5.1 represents a grain indenting in workpiece material and producing a 
groove of depth t . In the absence of friction at the surface between the spherical grain and 
the workpiece, the magnitude of the force required to indent the workpiece is constant and 
independent of the direction in which it is loaded [41,54], Therefore, radius b of the 
projected area of indentation is given by. 


(5.1) 
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where. H, is hardness of workpiece material and h is the radius of the projected area of 
indentation 



Fig. 5.1 Schematic diagram of indenting spherical abrasive grain 


sin^' 


From geometry of Fig 5 2 


2b _2jt(dg-t) 



cos^ 




(5 2) 
(5.3) 


where, 6 defines the line of action of the indentating force. Depth of indentation ‘ t ’ has 
been derived in chapter 3 and is given by following equation. 




1 

kv 


2 i 




(5.4) 


Tangential component of the chip formating force 
terms of force as given below 

F,g = F^g sin 6 


Ffg (Fig.5.2) can be written in 

(5 5) 


Substituting values of F^g and sin^ from eqns (5.1) and (5.2) respectively in eqn. (5.5), 

2^Hdpr) 

K V 


F,. 


tg 
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= — - (using eqn 3.59) 

3 3 

= (27C (5.6) 



Fig. 5.2 Abrasive grain producing a chip 

In practice, there will be friction along surface DC (Fig. 5.2). Assuming the mean 
coefficient of friction between grain and workpiece material to be / , the friction force will 
be /F,,. 

Hence, tangential component of friction force per active grit, = f F^g cosO 

= f f 1 -^- — - (using eqn. 3.59) 

^g 

= fn H^tid ( 5 . 7 ) 

Total tangential force per active grain (F,^ ) is equal to the sum of chip formation 
force and frictional force per active grain. Hence, 

3 3 

= (271- + fx H^tid.-D^diZFl 

d„ 


(5.8) 
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The specific energy in machining and abrasive processes is defined as the energy 
expanded per unit volume of material removal 

Specific energy, u = Vf /{A' Vf ) (5 9 ) 


where, Vf is velocity of flow of media along workpiece surface. It is also assumed to be 
equal to chip flow velocity. The area A' of groove generated, which corresponds to the 
cross hatched area in Fig 5 1, is given by the following equation 


d„ 


A' = . ,/7^( 

4 dg ^ ^ 2 


-t) (fromeqn 3.62) 


Total contact area between workpiece and media = n d^l^ 


(5 10) 
(5.11) 


where, d^ )s diameter of cylindrical workpiece and is length of workpiece 

Let n is the number of all simultaneously acting grits per unit contact area (active 
grain density) of media and workpiece. Then, total no. of grains simultaneously acting on 
workpiece surface = k nd^l^ (5 12) 


The tangential force acting between the workpiece surface and media is equal to 
the total of the tangential forces of all active grains within the contact area between the 
media and workpiece. 

Therefore, total tangential force = tangential component of cutting force + tangential 
component of friction force 

3 3 OA 

/.; = (JT « rf, /,) [(In H,) (rfj -lyid^ + fTt H^t (rfj - 0 ‘ ] 


“a 

In the above eqn. (5 13), dg, d^, and are the parameters that can be 


determined in advance. 

Total number of grains passed per stroke =2 tv (5.14) 

where, and areinner radii ofmedia cylinder and workpiece respectively, and I, is 


stroke length. 

As it is assumed that the path traced by an individual grain is straight line, volume 
of material removed by a grain is equal to cross sectional area of groove produced 
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multiplied by length of contact along the workpiece surface To simplify the analysis, 
length of contact is assumed to be equal to the length of workpiece surface. Therefore, 

Volume of material removed per unit time = (2 tc ry^nl^-^){A' 1^^) I {IJ (5 15) 

' w 

where Vp is the piston velocity. 

Energy required per unit time = specific energy x volume of material removed per unit 
time 

= u {2 TV r^nl,^)(A'lJI(l,IVp) (5.16) 

T 

'w 

From continuity equation, 


Vf _ 


V 


P 



(5.17) 


For the situation, where energy is dispersed uniformly over the workpiece contact 
area (2 n 1^), the thermal flux [energy (time)'*(area)'^] will be 

q==[u {2 7 c r^nl^\){AlJI{lJvp)]l{2nr^l^) 

= [unA' Vp -^] (5.18) 

For AFM, the specific energy is generally much larger than for other metal cutting 
operations. In practice, the specific energy can be used for the calculation of surface 
temperatures. From a more fundamental view point, an understanding of the origin of 
AFM energy can provide the basis for physical description of the mechanics of the 
process. 


5.3 Thermal Analysis of AFM Process 

The thermal behavior of workpiece is analyzed by one dimensional unsteady state 
heat transfer theory. Following assumptions are made to simplify the thermal analysis of 
AFM process' 

1 . The quantity of the heat flowing into the workpiece and media is constant during a 
machining stroke. 
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2 The heat flows along the thickness direction of workpiece and media. The unsteady 
temperature distribution in the thin workpiece, which is heated on one side by 
machining action can be approximated as one dimensional. Conduction in the 
direction of motion is neglected. The value of Peclet number obtained is of the order of 
1 O'* for all the cases considered. Hence, the conduction in the direction of motion can 
be neglected [55] 

3 The composite properties of media are considered in the evaluation of temperature 
distribution 

4. Thermal properties of media and workpiece are constant with time and space. 

5 Media is cooled to initial temperature in media cylinder after each stroke. 

The energy expanded in AFM is almost entirely dissipated as heat in the media- 

workpiece contact area, and causes the abrasive media and wokpiece temperature to rise. 

The heat generated in the AFM zone (Fig. 5.3) is removed from the region in the variety of 

ways 

1 It may conduct into the workpiece, where it is finally convected away by surrounding 
air. 

2 It may conduct into the media, where it is carried away by the motion of media. 



Fig. 5.3 Model of abrasive flow machining zone 
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Workpiece is assumed as one dimensional semi-infinite plate for the heat transfer 
analysis. For the same uniform initial temperature distribution, the surface of the 
workpiece is suddenly exposed to a constant heat flux . 

For constant properties of workpiece, the differential equation for the temperature 
(T) distribution in workpiece as a function of time T (z, r) is given by the following 

equation: 

^ = _L^ (5.19) 

dz^ dr 

where is the thermal difflisivity of the workpiece material and r is time. 

The initial and boundary conditions are 
T(z,Q) = T, 

^ Z ^ ^ 

7 ( 00 , r)=T, (5.20) 

where 7^ is the ambient temperature and the is thermal conductivity of workpiece 

material. 

This problem has been solved by the Laplace transform technique [56]. 

The solution for this case is 


T-T 



Aa^r 


(5.21) 


where is the heat flow to the workpiece per unit time and per unit area, z is the 
distance from the surface along the thickness, and is the thermal conductivity of 

workpiece material. 

The Gauss error function is defined as 

y. 2 x/2(VSt) 2 

erf — 7= = -Y= J ^ dq (5 22) 

2-yJar -sjTr o 

The energy equation for the moving media is considered at steady state to avoid the 
complication of the numerical part. It is also reasonable since in AFM time required for 
moving media to complete one stroke is very less in comparison to total time required for 
machining a workpiece. As workpiece is fixed, the temperature distribution in the 
workpiece is analyzed with respect to time. 
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The energy equation for media can be witten as 


(P^p)m — = k.. 


(5.23) 


where, Vm is the velocity of media, pj^, k^, Cp^ are the density, thermal conductivity and 

specific heat of media, respectively. 

The boundary conditions for T{x,z) are 

no,z) = T, 

P T , 


- k„ 


lz=0 


d z 

Tix,^) = T, (5.24) 

The solution to the differential equation (5.23) and boundary conditions, in non- 
dimensional form are 

r2 


T=2 ((— )2 exp (^) - Z erfc (-^) 


(5.25) 


T is non-dimensional temperature, and it is given by 


r = 


iT-T,) 


{qmKpCp'^)m] 

Other non-dimensional parameters are given by 


z v„ 




xv„ 




«c 

The dimensional form of eqn (5.25) is 


T-T 


[2 (^^)2exp 


^ 2 ^ V ^ 


(iZsL) erfc 


f(— 

2 xa^ 


V 


(5.26) 


(5.27) 


Let T be surface temperature at T(x,0) . From the eqn (5.27), the surface temperature is 


evaluated and is given by following expression. 

1 
N 

X 


r 


T,-T,=2q„ 


nik pc v\ 


(5.28) 


P vmy 


The average surface temperature on the workpiece surface can be expressed as follow. 
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(5.29) 


The surface temperature of workpiece obtained by the eqn. (5.21) can be written as follow 

(5.30) 


T,-T, 


_ 2^^ \a^ r 




For the common average surface temperature (equating eqns (5.29) and (5.30)), 

9/, V _ 2g>, r 

\n:(kpCpV)^] k^ i n 

Let R be the fraction of the total heat flux q entering into the workpiece, that is 
generated along the workpiece and media interface. 

^ = + =R^ + (1-R)9 (5.32) 

Then, after simplification of the eqn. (5 32) the fraction of heat entering the 
workpeice can be obtained by the following expression. 

1 


R = 


L 


''a 


ipc.kv) 


mj 


(5.33) 
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5.4 Evaluation of Properties of Media 

In order to evaluate temperature distribution, the properties of media must be 
known. A simple weighted mean model for specific heat capacity based on the relative 
properties (by mass) of each constituent in media, is used to determine the variation of this 
property in the mixture. 

= (1- C/100) + (C/ 100) (5.34) 

where, C is the percentage concentration of abrasives in the media. c_, and c are the 

P pQ 

specific heat of carrier (polyborosilixane) and abrasives respectively. 

The density of media p^ in terms of abrasive density p^ and carrier density p^ is 
given by the relation. 
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(5.35) 


[C>,+(100-C)pJ 

The thermal conductivity of media [22] can be obtained by the following 
expression 


1 




(5.36) 


S] = 0.8061 d 


S2 = Sj 


OOP -Q Pa 
C Pa 


+ 1 


S3 =S2- Si 

where, and are thermal conductivity of carrier and abrasives, respectively. 


Flow of media 



Fig. 5.4 Experimental arrangement for temperature measurement 


5.5 Experimentation 

Experiments have been conducted on the AFM set-up designed and fabricated in 
house Internal surfaces of mild steel workpieces have been finished using media 
(polyborosiloxane carrier mixed with SiC abrasives). The experimental arrangement for 
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measurement of temperature is shown in Fig. 5.4. Temperature at the workpiece surface 
during AFM has been measured using copper-constantan thermocouples. A multi-channel 
digital temperature indicator having resolution of 0.1 was used for recording the 
temperature. The conditions, under which AFM experiments have been carried out, are 
given in Appendix C Details of experiments are discussed in chapter 2 


5.6 Results and Discussion 

The conditions of machining and values of various parameters, given in Appendix 
C, have been considered for the evaluation of forces, specific energy and heat flux The 
values of active grain density n , reported in chapter 4, have been used in the calculations. 
As machining action in AFM compares to grinding, the coefficient of friction between 
workpiece surface and abrasives is assumed to be 0.2 This is also cited in the literature 
[41] for calculation of specific energy in grinding. It is obvious that specific energy is a 
function of the average cross sectional area of a cut. The average cross sectional area of a 
cut seems to be smaller in AFM than grinding With work material of low hardness the 
cross sectional area of cut is larger with the same normal force. The specific energy in 
AFM (Fig. 5.5) is near to the range of specific energy in grinding (5-60 J/mm^) [41] for 
mild steel workpiece. Generally it is higher in AFM due to lower depth of indentation. Fig. 
5.6 indicates that specific energy decreases with increase in normal pressure due to higher 
depth of indentation. The normal pressure of indentation have been determined from the 
extrusion pressure applied in AFM, as explained in chapter 3 It is interesting to note that 
specific energy remains almost constant with change in abrasive mesh size (Fig. 5.7). The 
similar results have also been reported in case of grinding by Kannappan and Malkin [57]. 

A fraction of generated total heat initially going to the workpiece is 96.5% Its 
value reduces with the time of machining (Fig. 5.8). As a major portion of the heat enters 
into the workpiece, rise in temperature of the workpiece is much higher than rise in 
temperature of the media. It is due to much higher thermal conductivity of workpiece 
{k^=60.5 W/m-K) as compared to the media. The thermal conductivity of media is 0.38 

W/m-K for 66 % concentration of abrasives with average mesh size of 55. Temperature of 
workpiece increases with number of cycles (Fig. 5.9). Initially rate of rise in temperature is 
more, and it decreases with number of cycles. Measured values of rise in temperature are 
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lower than calculated ones due to simplification of analysis to one dimensional heat flow 
It is also observed (Fig 5 10) that rise in temperature increases with extrusion pressure, as 
higher energy at higher pressure results into increase in heat flux at workpiece surface 



Fig. 5.5 Variation of specific energy with hardness of workpiece material. 



Normal pressure, MPa 


Fig. 5.6 Variation of specific energy with normal pressure of indentation 
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Fig. 5.7 Variation of specific energy with mesh size of abrasive grains 
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Fig. 5.8 Variation of percentage of heat entering the workpiece with machining time 
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Fig, 5.9 Change in temperature of workpiece with number of cycles 
(C = 66 %, = 50 bar, = 0.905) 



Fig. 5.10 Change in temperature of workpiece surface with extrusion pressure 

( C = 66 %, N = IS, Re = 0.905) 
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5.7 Conclusions 

The present analysis for specific energy and tangential force is a step forward in 
understanding the mechanism of AFM process. The specific energy for AFM (10-110 
J/mm^) process is partially overlapping the range of specific energy for grinding (5-60 
J/mm^). This justifies the inference of researchers that machining action in AFM compares 
to grinding It is interesting to note that specific energy remains constant with abrasive 
mesh size, but its value is higher for higher hardness of workpiece material 

A simple model has been developed for the heat transfer in AFM process 
ccusidering the heat flow to the workpiece and media A fraction of the heat entering the 
workpiece and media can be determined. The model predicted the workpiece temperature 
with reasonable accuracy However, The present analysis can be improved upon to give 
more realistic temperature distribution, by modeling the problem as a three dimensional 
unsteady state heat transfer problem 



Chapter -6 

Application of Neural Networks and Genetic Algorithms 
to AFM Process 


6.1 Neural Networks in Manufacturing 

The success of automated manufacturing relies to a large extent on the development 
of computer-based learning schemes that are able to code operational knowledge. The 
complex requirements of many modem applications demand for flexible/adaptable devices 
that can quickly react to changing circumstances as solution. Today's complicated machine 
tools, manufacturing cells and systems require systematic methods of monitoring, control 
and fault diagnostics. Neural network systems promise this capability. Process modeling 
and optimization are very important issues in manufacturing engineering. Abrasive flow 
machining process is too complicated to warrant appropriate analytical models and most of 
the time, analytical models are developed based on many simplified assumptions, which 
contradict reality. More importantly, it is sometimes difficult to adjust the parameters of the 
above mentioned models according to the actual situation of the machining process [58]. 
Because of the complexity of AFM process, the optimization as well as the optimal control 
is difficult to perform. Neural networks provide significant advantages in solving machining 
problems that require real-time encoding and interpretation of relationships among variables 
of high dimensional space. Neural networks have been shown to be effective for modeling 
complex non-linear processes. Unlike techniques such as non-linear regression, neural 
networks do not require a priori assumption of the functional form of the model (linear, 
polynomial and so on). Therefore, neural networks which can map the input/output 
relationships and possess massive parallel computing capability, have attracted attention m 


research. 



120 


Application of Neural Networks and Genetic Algorithms to AFM 


Recently, neural networks have been applied in precision manufacturing to predict 
some machining phenomenon as well as machining quality. Chiyssolouris and Guillot [59] 
modeled the machining processes by a multiple regression method and neural networks, and 
concluded that neural networks are superior to conventional multiple regression methods 
Rangwala and Domfeld [60] presented a scheme that used a multilayered perceptron neural 
network to model the turning process and an augmented Lagrange multiplier method to 
optimize the material removal rate (MRR). A simple neural network model for grinding 
mode identification and surface quality prediction in grinding of silicon nitride is also 
established by Zouaghi and Ichida [61], Sathyanarayan et al [62] used a neural network 
model to study creep feed grinding of super alloys, but the optimization was done 
analytically using an off -line multi-objective programming technique Liao and Chen [58] 
showed how multi-layer perceptron (MLP) networks can be used to model and optimize 
grinding processes. In their work, a generalized four-layer perceptron network was trained 
with the standard back propagation (BP) algorithm to model the process and a Boltzman 
factor was integrated with the BP algorithm for process optimization 

It is the intent of this chapter to realize ways to evaluate accomplishment and 
performance of neural network in modeling and optimization of AFM process. In this 
chapter, architecture of neural network for prediction of material removal rate and surface 
roughness in AFM process is presented. First, a generalized back-propagation neural 
network is used to establish the process. model. The modeling phase is followed by an 
optimization phase, during which the network predicts the input conditions to be used by 
AFM process to maximize material removal rate subject to appropriate operating 
constraints on surface roughness and process parameters. The optimization of the AFM 
process is also done using genetic algorithm. 

6.2 Overview of Neural Networks 

Neural networks are developed to model the way in which the human brain 
performs a particular task, or process information. A neural network is a massively parallel 
distributed processor that has a natural propensity for storing experiential knowledge and 
making it available for use. The central motivation underlying the development of artificial 
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acquired and stored over time through the use of adaptive learning algorithms [63], It 
resembles the brain in two respects : 

1 Knowledge is acquired by the network through a learning process. 

2 Interneuron connection strengths known as synaptic weights are used to store the 
knowledge. 

A neural network is defined by three features, topology, functionality, and learning 
Topology refers to the number of layers, number of nodes in each layer, and the way nodes 
are connected. Functionality refers to the transfer function and discriminatory function (if 
any) of each node, and the cost function of the network outputs Learning refers to the 
learning algorithm and the values of learning parameters (e.g. learning rates, and 
momentum). Learning rules specify an initial set of weights and indicate how weights 
should be adapted during use to improve performance of the neural network system 

Neural networks can be classified differently according to different criteria : the 
type of learning (supervised, unsupervised, self supervised), the mode of operation (on-line, 
off-line), the network interconnection architecture (feedforward, feedback, recurrent), the 
number of layers, the learning algorithm (Hebbbian, back-propagation, resonant, Widrow- 
Hoflf, simulated annealing), and the type of data processed (binary, discrete, continuous, 
fuzzy). 

The most popular neural network topology is feedforward networks, viz. multilayer 
perceptron (MLP). The back-propagation (BP) network is the most commonly used neural 
network because there exists a mathematically strict learning scheme to train the network 
and guarantee mapping between inputs and outputs. Therefore, MLP neural network with 
back-propagation algorithm is used in the present study for modeling and optimization of 
AFM process. 


6.3 Back-Propagation Neural Network Model for AFM Process 

Recently, it has been shown that multilayer perceptrons (MLP) neural networks are 
suitable for the modeling of manufacturing processes [58]. Multilayer perceptrons have 
been applied successfully to solve some difficult and diverse problems by training them in a 
supervised manner with a highly popular algorithm known as the error back-propagation 
algorithm. A multilayer perceptron has three distinctive characteristics [64]. 
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1 The model of each neuron in the network includes a non-linearity at output end. The 
important point to emphasize here is that the non-linearity is smooth (i e , differentiable 
ever 5 where). 

2 The network contains one or more layers of hidden neurons. These hidden neurons 
enable the network to learn complex tasks by extracting progressively more meaningful 
features from the input patterns. 

3 The network exhibits a high degree of connectivity determined by the synapses of the 
network. 



Fig. 6.1 Structure of three layered neural network 

Back-propagation neural network [64] is usually referred to as feed forwarded, 
multilayered network with number of hidden layers (Fig. 6.1) trained with a gradient 
descent technique. This algorithm is based on the error correction learning rule. Basically, 
the error back-propagation process consists of two passes through the different layers of 
the network • a forward pass and a backward pass. In the forward pass, an activity pattern 
(input vector) is applied to the sensory nodes of the network, and its effect propagates 
through the network layer by layer Finally, a set of outputs is produced as the actual 
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response of the network During the backward pass, all synaptic weights are adjusted in 
accordance with the error correction rule Specifically, the actual response of the network 
is subtracted from the desired (target) response to produce an error signal. The synaptic 
weights are adjusted so as to make the actual response of the network move closer to the 
desired network 

Let, i and j = Different neurons of the network 

w^ , jt = Synaptic weight vector between neuron yin layer (A:-l) to neuron i in 

layer k . 

Vjk = Net internal activity levels of neuron j in layer k . 
yj i^ = Function signals of the neuron j in layer k . 

The internal activity level ^ is 


''j.k 




( 6 . 1 ) 


;=1 


where, p is the total number of inputs applied to neuron j . 

Assuming use of a logistic function for the sigmoidal nonlinearity, the function 
(output) signal of neuron j in layer k is 

] ( 6 . 2 ) 


yj,k= 


[l + exp(-v^.,*)] 


If neuron j is in the first layer, yj,i=Xj 
and if neuron J is in the output layer, yj^k = 

where Xj is the y th element of the input vector x and is a function signal of the 


neuron j in the output layer. 

The average of the instantaneous sum of the squared errors of the network 
1 ± 


ze; 

2 ■' 




■o,r 


( 6 . 3 ) 


where q is the total number of output nodes, rf, is the ;th element of the desired 
response, d . 

The back-propagation algorithm applies a correction Aw,,, to synaptic weights 


w, , , as given below. 
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where 77 is a constant and it is called learning rate parameter of the back propagation 
algorithm. 

The local gradient 5 of the network, obtained by proceeding backward layer by 
layer, is given as 

(6.5) 

for neuron j in the output layer. 

5j,k = yj,S- yj,k)'ZS ]M\'^ j,>M\ 

1 

for neuron j in the hidden layer k . 

A simple method of increasing the rate of learning and yet avoiding the danger of 
instability is to modify the delta rule by including a momentum term [65]. Then, 

Aw,_^ (n) = a w,_^ (n - 1) + 77 S^y^ (n) (6.7) 

where a is usually a positive number called ‘momentum constant’, n is the number of 
training pattern presented to the network. 

The computation is iterated by presenting the new set of training examples to the 
network until the jfree parameters of the network stabilize their values and the average 
squared error computed over the entire training set is at a minimum or acceptably small 
value. 

6.4 Optimization Using Neural Network 

Nonlinear constrained programming is a basic tool in the systems where a set of 
design parameters is optimized subject to inequality constraints. The goal of this section is 
to propose a parallel structure, one that parallels the conventional optimization algorithm to 
solve this problem The key point in the optimization phase is to find the partial derivatives 
of the cost function with respect to the input data, not the partial derivatives of error with 
respect to the parameters of the network as in the modeling phase. Further, the process 
optimization is assured by minimizing a cost function L subject to the given constraints. 

The constraints on input variables are 

x,(min) < X, ^ X, (max) (f = 1 ,/?) 


( 6 . 8 ) 



Application of neural networks and Genetic Algorithms to AFM 


125 


and the constraints on the output variables are 

gj(x) = (cij-yj )<0 (J = l,q) (6.9) 

where x, (mm) and x,(max) are the minimum and maximum values of the input variables 
X, , respectively In AFM process, the constraint is on surface roughness. 

In the simulation model presented below, the augmented Lagrange multiplier 
(ALM) method [66] is adapted to solve the above problem. 

According to ALM method, the pseudo-objective function is given as 

A( X, 1 , tp) = L(x) + ^ {gj (x) + sj} + Tp {gj (x) + sj}^] (6. 10) 

where s^ are slack variables added to inequality constraints; A, is the Lagrange multiplier; 
and rp is a penalty parameter. 

A( X, A , rp) = L(x) + /y + fp ^f] (6. 1 1) 

j ■' 

. 

where, = max { ^/x) , 

^ r_ 


The upgrade formula for is 


( 6 . 12 ) 

The differentiation of the pseudo-objective function (given above) is obtained as 
5 A 


dx, 


5L o j. 


( 6 . 13 ) 


where Ry = 1 


if gjix) >- 



= 0 otherwise 

The partial derivatives of L and g can be calculated by using back-propagation 
algorithm. The amount of change in initial input values is proportional to the derivative of 
the objective function 


dA 




ayj 


( 6 . 14 ) 


where is a parameter, which decides the rate of change in initial input values. The 
gradient o-,,* of the network, obtained by proceeding backward layer by layer, is given as 
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CT,, = yj,, (1 - )] for ( ^ - 1 )th layer (6. 1 5) 

CTa-i = IK, w,, .,_i(l-;/,.fc_,) for lower layers. (6.16) 

For the first layer, (6-17) 

dy, ax, 

The backward pass starts at the output layer by passing the change in objective 
function signals leftward through the network, layer by layer, and the recursively computing 
gradient or, for each neuron. The recursive process permits the input of the network to 
undergo change in accordance with the analysis. The process is continued, until the input of 
the network stabilizes their values and the change in objective function value is at a 
minimum or acceptably small value. 

In AFM process, the cost function L(x), which is to be minimized, is negative of 
the material removal rate and constraint g(x) is on the surface roughness, which should be 
less than or equal to maximum allowable surface roughness. 

The solution to the above problem is a locally optimal set of process input 
variables. The optimization by neural network is simply an optimization procedure with two 
important differences: no assumption of an input-output model is made, and the 
optinaization procedure is made in parallel. 


6.5 Optimization Using Genetic Algorithms 

Optimization of AFM process has also been carried out by genetic algorithm to 
compare the results of optimization by neural network. Even though the traditional or 
classical optimization algorithms [66,67] have been extensively used in many engineering 
design and decision-making problems, there are a number of shortcomings. Most of the 
popular methods require the gradient information, which may not be easy to calculate (or 
not available at all) in many real-world problems. Moreover, most of the traditional 
methods may converge to a locally optimal solution which may be very different from the 
actual global optimal solution. Further, none of these methods can be applied to a wide 
variety of problems. In order to alleviate some of these problems, a couple of non- 
traditional search and optimization algorithms - genetic algorithms and simulated annealing 
- are finding wide spread applicability in engineering design and decision making problems. 



Application of neural networks and Genetic Algorithms to AFM 


127 


These methods work according to the principles of natural phenomenon and found to have 
successfully solved many non-linear programming (NLP) problems to a global optimality 
[68] In the next section, we describe the introduction and working principles of genetic 
algorithms 

6.5.1 Introduction to genetic algorithms 

Genetic Algorithm (GA) can be understood as an “intelligent” probabilistic search 
algorithm which can be applied to a variety of combinatorial optimization problems. The 
theoretical foundations of GAs were originally developed by Holland [69], GAs are based 
on the evolutionary process of biological organisms in nature [70]. During the course of 
evolution, natural populations evolve according to the principles of natural selection and 
“survival of the fittest” Individuals which are more successful in adapting to their 
environment will have better chance of surviving and reproducing, whilst individuals which 
are less fit will be eliminated. This means that the genes from the highly fit individuals will 
spread to an increasing number of individuals in each successive generation. The 
combination of good characteristics from highly adapted ancestors may produce even more 
fit offspring. In this way, species evolve to become more and more well adapted to their 
environment. 

6.5.2 Genetic algorithms - working principle 

As described in the preceeding section, Genetic Algorithms are search and 
optimization procedures that are motivated by the principles of natural genetics and natural 
selection. Some fundamental ideas of genetics are borrowed and used artificially to 
construct search algorithms that are robust and require minimal problem information. 
Consider a typical unconstrained single variable optimization problem, which can be stated 

as [70]: 

Maximize: fix) (6-18) 

Variable bounds: <x^ 

In order to use GAs to solve the above problem, the variable x is typically encoded 
into a string or chromosome structure which represents a possible solution to the given 
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problem. Binary coding is mostly used and the length of the string is usually determined by 
the accuracy of the solution desired For example, if five bit binary strings are used to code 
the variable x then the string (0 0 0 0 0) is decoded to the value Xmin, the string (1 1 1 1 l)is 
decoded to the value Xnmx, and any other string is decoded to a value in the range (Xm,n, Xmax) 
uniquely. Thus, with five bit strings used to code the variable x, there will be a total of 2^ or 
32 different strings that are possible, and the accuracy between two consecutive strings will 
be (x„,ax - Xnu„)/32. If more accuracy is desired longer strings can be used. With a known 
coding, any string can be decoded to a x value, which can be used to find the objective 
function value. A string’s objective function value (/(x)) is called as the string’s fitness. 
A pseudo code for a simple genetic algorithm is given in Figure 6 2. GAs begin with a 
population of strings (individuals) created at random [70]. Fitness of each individual is 
evaluated with respect to the given objective function. Then this initial population is 
operated by three main operators - reproduction, crossover, and mutation - to create 
hopefully a better population Highly fit individuals or solutions are given opportunities to 
reproduce by exchanging pieces of their genetic information, in a crossover procedure, with 
other highly fit individuals. This produces a new “offspring” solutions (i.e. children), which 
share some characteristics taken from both the parents. Mutation is often applied after 
crossover by altering some genes (i.e. bits) in the strings. The offspring can either replace 
the whole population (generational approach) or replace less fit individuals (steady-state 
approach). This new population is further evaluated and tested for some termination 
criteria. If the termination criteria are not met, this reproduction-crossover-mutation- 
evaluation cycle is repeated until the termination criteria are met. One cycle of these three 
operators and evaluation procedure is called a generation in GA terminology. Before 
discussing the application of genetic algorithms to the present work, let us briefly discuss 
these three GA operators for a better understanding of the GAs. 

Reproduction selects the good strings in a population and forms a mating pool for 
the subsequent crossover. There exist a number of reproduction operators in GA literature, 
but the essential idea is that above average individuals are picked from the current 
population and duplicates of them are inserted in the mating pool. The commonly used 

reproduction operators are proportionate selection with roulette-wheel and tournament 
selection. 
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Begin 

Initialize population; 

Evaluate population; 
repeat 

Reproduction; 

Crossover; 

Mutation; 

Evaluate population; 
until (termination criteria); 

end. 

Fig. 6.2 A Pseudo-code for a Simple Genetic Algorithm (SGA) 

Crossover operator is applied next to the strings in the mating pool. Here, two 
strings are picked at random from the mating pool and some portion of the information is 
exchanged between the strings. For example, in a single point crossover operator, this is 
performed by randomly choosing a crossing site along the string and by exchangmg aU bits 
on the right side of the cross site. In order to preserve some good strings found in the 
mating pool, not all striogs in the population are used in crossover. If a crossover 
probability of p. is used, then (paXlOO)% of strings are used for crossover and the rest of 

the strings are simply copied to the new population. 

Mutation changes a 1 to a 0 and vice-versa with a small mutation probability p„ 

The need for mutation is to keep diversity in the population. For example, if in a particular 
position along the string length all strings in the population have a value 0, and a 1 ts 
needed in that position to obtain the optimum then neither reproduction nor crossover w, 
be able to create this. But mutation introduces some probability of turning that 0 mto a 1. 
Furthermore, mutation may be found useful for local improvement of a soluuon. 
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6.5.3 Optimization of AFM process 

The problem of optimization of AFM process can be described as minimizing or 
maximizing an objective function subject to a certain set of constraints. In the present 
analysis, the objective function is maximization of material removal rate (MRR), which is 
defined in terms of process variables. The process variables in the present case are average 
piston velocity (Vp), percentage concentration of abrasives (C), mesh size of grains (A/) 
and number of cycle (N) In order to limit the search space, constraints in the form of 
lower and upper bounds on the process variables are imposed. 

The equations given by Adsul [19], for predictions of MRR and surface roughness, 
were used here for optimization of AFM process. These models have been developed by 
multi variable regression analysis (MVRA). To estimate the regression coefficients, the 
method of least squares has been employed. The material removal rate (MRR) in g/min can 
be expressed as [19] 

MRR= 5.285 x lO''® (Vp)^®"®^ (C)'-"’® (Nj" (6.19) 

The surface roughness (i?^ ) in jum can be expressed as 
R, = 2.8275 X 10^ (Vp)'‘®“^ (C)'^-^“^ 20) 

The problem of optimization of AFM process can be described as maximizing MRR 
subject to a certain set of constraints on surface roughness and input variables. 

In order to use GA the constrained optimization problem is stated as follows; 

Maximize: MRR (6.21) 

Subject to: (6.22) 

where and x" are the lower and upper bounds on the process variable x, . 

A Simple Genetic Algorithm code was used in the present study. The steps 
involved in the optimization using GA of AFM process can be stated as follows: 

STEP 1: The GA parameters are initialized. This involves specifying the population size, 
maximum number of generations, string length of each variable, crossover and mutation 
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probabilities etc Upper and lower limits on each of the process variables are also specified 
The generation number is initialized to 0 

STEP 2; An initial feasible random population is generated 

STEP 3 . The fitness of each individual in the population is evaluated 

STEP 4 . Once the fitness of all the individuals is available, GA operations are performed, 

e g. reproduction, crossover, and mutation. Now, a new set of process variables is created 

which possibly is better than that of the previous generation. This completes one 

generation. 

STEP 5: The generation number is incremented. If the current generation number is 
greater than the maximum number of generations, then the process is terminated. 
Otherwise, process is repeated from step 3. 

\ 

6.6 Experimentation 

Machining experiments were carried out by Adsul [19] on a set-up of AFM machine 
using aluminum as workpiece material. The media was a mixture of putty and silicon 
carbide abrasive particles. Experiments were designed and conducted for various sets of 
process variables, i.e. average piston velocity, abrasive concentration, abrasive mesh size 
and number of cycles. For performing these experiment compositions of different abrasive 
mesh sizes (100 to 240) and having different abrasive concentrations (33 to 45 %) were 
prepared. Percentage concentration of abrasives is defined as the ratio of weight of 
abrasives and total weight of abrasive media (abrasives +carrier) multiplied by 100. The 
piston velocity range chosen was from 40 to 85 cm/min. Average material removal rate 
(MRR) over test duration was measured. In these experiments, the effect of various process 
parameters on surface roughness and material removal rate were studied. The data 
collected from these experiments have been used to train the neural network and to test the 
validity of the present model. The details of experimental results used in this paper are 
given in Appendix D. 
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6.7 Simulation Procedure 

The objective of the simulation using neural network was to first have the system 
learn the appropriate mappings between input and output variables by observing the 
training samples. The trained system was then used to determine the input conditions that 
maximizes material removal rate (MRR) subject to certain constraints. The three layer 
back-propagation with four inputs, two outputs and nine hidden nodes was employed for 
neural network. The network was trained using 25 samples that span the allowable ranges 
of input variables. Due to limited number of data available in the work by Adsul on AFM 
experiments restricted us to take only 25 data points in this study. The inputs were 
average piston velocity (Vp ), percentage concentration of the abrasives (C), abrasive mesh 

size {M ) and number of cycles {N), and outputs were material removal rate (MRR) and 
surface roughness A random generator was used to initialize the values of the 

learning parameter. In order to decide the structure of neural network, the rate of error 
convergence was checked by changing the number of hidden layers and also by adjusting 
the learning rate and momentum rate. As a result, a neural network with nine neurons in the 
hidden layer was adopted for storing knowledge in the form of weights between neurons. 
The simulation results for the initial training indicate (Fig. 6.3) that it takes about 4000 
iterations for the root mean error to reduce to 2%. 

A crucial problem in back-propagation is its generalization ability. There is no 
certainty that the network successfully trained on the given samples provides desired 
input/output associations for untrained pattern as well. Concerning this problem, some 
authors [71] showed experimentally and mathematically that the generalization capability 
could remarkably be enhanced by training the network with noise injected inputs. The 
addition of random noise allows the system to jump out of the local minima that may trap it 
in a high error state. In the present simulation, the network was trained by both noiseless 
data and noise injected data. Fig. 6.4 shows five different cases studied, and explanation of 
different abbreviations used. For example, abbreviation NioNi indicates that network was 
trained with noisy data (le noise is added in both inputs and outputs), and prediction is also 
for noisy input. In the Fig. 6.4, information after first capital letter follows for training data 
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and that after second capital letter follows for prediction data It is seen that the root mean 
percentage error is almost the same for various cases for low percentage of noise 
However, it increases rapidly with the increase in percentage noise except for the case of 
NiNi which means training with noisy input variables and prediction for noisy input It is 
observed that for the same parameters of the network, the total error with 2.5 percentage 
random noise injected inputs is lesser than the error obtained with noiseless inputs It is also 
found that training with 2 5 % noise injected inputs provides enhanced learning compared 
to noiseless input (Fig. 6.4) Variations is observed in the response of the network during 
learning but the difference between the two is not significant after about 4000 iterations. 
The number of iterations required to achieve same amount of error in noise injected 
learning is less than the number of iterations required in noiseless learning. For example, for 
2.0 % error, number of iterations required with noiseless input are 3710, while with 2.5% 
noise injected input, it is only 3390 (Fig 6.4). Thus, noise control method consistently 
produced faster learning than standard back-propagation 



Fig. 6.3 Variation of root mean percentage error with number of iterations 
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Percentage Noise 


Fig. 6.4 EfTect of noise in back propagation learning. 

NioNi ■ Training with Noisy input and output, and prediction using Noisy input. 
NioNl : Training with Noisy input and output, and prediction using Noiseless input 
NiNi • Training with Noisy input, and prediction using Noisy input. 

NiNl : Training with Noisy input and prediction using Noiseless input. 

NlNi . Training with Noiseless input and prediction using Noisy input. 

6.8 Results and Discussion 

The experimental results have been compared with the results computed using multi 
variable regression analysis (MVRA) (eqns. (6.19) and (6.20)), and those obtained by 
neural network discussed above. The MRR and surface roughness results predicted by 
simulation using neural network show a good agreement with the experimental results and 
the results obtained by the MVRA models (Fig. 6.5-6.11) for a wide range of operating 
conditions. Variation of surface roughness with number of cycles (Fig. 6.9) indicates that 
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error in prediction is higher at lower number of cycles At higher value of number of cycles, 
predictions by neural network are in good agreement with experimental values. A 
comparison of results predicted by the present neural network model and regression model 
also revealed that absolute percentage error in prediction of untrained data by neural 
network model was from 0.25 % to 8.95 %, whereas it was from 0.09 % to 25 % in results 
predicted by multi variable regression analysis 

Next, the network has been used to determine a set of optimal inputs. In the present 
case of optimization of AFM process, the cost function, L(x) which is to be minimized, is 
the negative of material removal rate (MRR) and constraint ^(x) is on the surface 
roughness obtained, which should be less than or equal to the maxim um allowable surface 
roughness, Ramax ■ The cost function is given as 


L(x) = - MRR 

subject to inequality constraints 

(6.23) 

g(x) R^^^ 

(6.24) 


According to the experimental set-up and workpiece used, the limits on the input 
variables are as follows : 


40.0 < Vp < 85.0 

33.0 < C < 45.0 

100.0 ^ M < 240.0 

20.0 < N < 120.0 

Table 6. 1 shows the optimal inputs synthesized by the network for different values 
of constraint on surface roughness. 

In order to optimize the present problem using GA, the following parameters has 
been specified by practice, to get optimal solutions with less computational effort. 
Population size = 50 
Max. number of generations = 200 
Total string length = 40 
Cross over probability =0.8 
Mutation probability = 0.01 
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The optimal inputs and corresponding values of outputs obtained by GA are 
presented in Table 6 2. The results obtained by neural network approach and genetic 
algorithm are very close. It indicates that the network possesses sufficient knowledge in 
order to perform optimization. 


Table 6.1 Optimization results of the network 


No 

D 

max 

jum 


optimal inputs 


optimal outputs 

Vp 

cm/min 

C 

M 

N 

K 

fi m 

MRR 

mg/min 

1 

07 

85.00 

44.92 

101.90 

20.78 

0.625 

0 716 

2. 

06 

84 60 

44 38 

100.44 

21.46 

0 590 

0 692 

3 

05 

83.59 

44 80 

101.17 

33.78 

0.500 

0.635 

4. 

04 

84 25 

44.56 

101.76 

54.40 

0.392 

0 538 

5 

0.3 

85.00 

44.88 

101.03 

119.70 

0.289 

0.516 


Table 6.2 Optimization results of the genetic algorithm 


D 

max 


optimal inputs 


optimal outputs 

jum 

''p 

C 

M 

N 

K 

MRR 


cm/min 




fi m 

mg/min 

0.7 

84 69 

44.45 

100.88 

21.09 

0.607 

0.697 

0.6 

84.82 

44.97 

101.02 

25.47 

0.553 

0 695 

0.5 

84.95 

45.00 

102.05 

33.00 

0.490 

0.669 

0.4 

85.00 

44.74 

100.43 

63.4 

0.370 

0.580 

0.3 

84.95 

44.96 

101.02 

111.13 

0.288 

0.525 


4. 

5. 




20 40 60 80 100 120 

Number of cycles 


Fig. 6.5 Effect of number of cycles on material removal rate 
(C = 45%, Vp = 51.5 cm/min, M= 100) 



30 34 38 42 46 

Percentage abrasive concentration 

Fig. 6.6 Effect of percentage abrasive concentration on material removal rate 
( Vp = 51.5 cm/min, M = 100, N= 80) 
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Fig. 6.7 Effect of abrasive mesh size on material removal rate 
(C = 45%, Vp = 51.5 cm/min, A = 50) 



Fig. 6.8 Effect of piston velocity on material removal rate 
(C = 45%, M = 100, N = 50) 
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Fig. 6.9 Effect of number of cycles on surface roughness 
(C = 45%, Vp = 51.5 cm/min, M = 100) 



Fig, 6,10 Effect of percentage concentration on surface roughness 
( Vp = 51.5 cm/min, M = 100, N = 80) 
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Fig. 6.11 Effect of piston velocity on surface roughness 
(G =45%, M= 100, A' = 50) 


6.9 Conclusions 

In this chapter, the effectiveness of using back-propagation neural network for 
process modeling and optimization of abrasive flow machining process has been 
demonstrated. Simulation results show a good agreement with experimental results for a 
wide range of machining conditions. The optimization results of the neural network also 
coincide well with the results obtained by genetic algorithm. Based on the simulation and 
optimization results, the following conclusions can be drawn . 

1 The possibility of using this neural network model for machined surface quality and 
material removal rate prediction for AFM process has been confirmed. 

2. It is observed that learning could remarkably be enhanced by training the network with 
noise injected inputs 
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3 Appropriately trained neural network successfully synthesizes optimal input conditions 
for AFM process. The optimal input conditions maximize the material removal rate, 
subject to appropriate process constraints, 

4 An important consideration is that process optimization can be performed in the 
absence of process models and purely by observations of experimental information. 

5 The discussed neural network system is fairly general and can be extended to other 
abrasive processes to improve machining efficiency. 



Chapter 7 

Conclusions and Scope for Future Work 


7.1 Conclusions 

The present thesis work comprises both experimental and theoretical 
investigations. Experiments have been conducted to study the influence of the process 
parameters on the process performance. A model has been proposed for analysing the 
media flow in AFM process through media cylinder and workpiece, and solved by finite 
element method. A model for the prediction of material removal and surface roughness in 
AFM has also been presented. The generation of surface profile and material removal in 
AFM have been simulated considering random distribution of abrasive grains in media and 
their interaction with workpiece surface profile during the process. Specific energy and 
thermal aspects of AFM have also been analysed. Further, the effectiveness of using back- 
propagation neural network for process modeling and optimization of AFM process has 
been demonstrated On the basis of above mentioned experimental and theoretical 
investigations, the following conclusions have been drawn. 

Experimental investigations 

. AFM experiments on mild steel workpiece using SiC abrasives mixed with 
polyborosiloxane carrier indicate that AFM successfirlly reduces the surface roughness 
of workpieces. AFM performance can be easily controUed by varying the process 

parameters. 

. An increase in extrusion pressure increases the depth of indentation of abrasives, and 
increase in percentage concentration of abrasives increases active grain density and 
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media viscosity, hence material removal and change in surface roughness are higher at 
higher extrusion pressure and concentration of abrasives in media 

• The slug length of media flow in a stroke increases with increase in reduction ratio. 
Hence, material removal and change in surface roughness increase with increase in 
reduction ratio. 

Analysis of material removal and surface roughness in AFM 

• A good agreement exists between the theoretical and experimental values of extrusion 
pressures at various values of piston velocities. It clearly indicates that the proposed 
flow model provides a reasonably good estimate of the flow characteristics involved in 
AFM process 

• At a particular extrusion pressure, piston velocity is higher at lower angle of tooling. 
But the length of tooling increases with decrease in angle of tooling, a . Hence, there is 
a need to determine the optimum angle of tooling for minimum extrusion power and 
machining time. 

• The piston velocity increases linearly with extrusion pressure for a specified reduction 
ratio. At higher value of reduction ratio, the rate of increase in extrusion pressure is 
also higher. The normal stresses on the workpiece surface increase with increase in 
reduction ratio Finite element analysis of media flow has been carried out to predict the 
normal stresses at the workpiece surface. 

• A model has been developed for prediction of material removal and surface roughness 
by applying abrasion theoiy. The model is based on accumulated plastic flow of 
workpiece material by repeated indentation by moving abrasive particles. Good 
agreement is found between the response surface analysis (RSA) and theoretical results. 

• Material removal increases with extrusion pressure, reduction ratio and percentage 
concentration of abrasives for a specified number of cycles. 

• Change in surface roughness (A/?<ar) value increases with increase in extrusion pressure, 

percentage concentration of abrasives, and reduction ratio for a specified number of 
cycles. 
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Simulation of surface generation in AFM 

• The generation of surface profile and material removal in AFM have been simulated 
considering random distribution of abrasive grains in the media and their interaction 
with workpiece profile during machining. The results of the proposed simulation and 
RS A are in good agreement, which justifies the use of this approach for predicting the 
expected surface characteristics in AFM 

• Two approaches to determine active grain density in media — first, optical microscope 
technique and second, stochastic simulation method have been proposed. The results of 
stochastic simulation method are in agreement with the data obtained by microscopic 
technique. 

• Active grain density increases with an increase in abrasive concentration, and for a 
particular concentration, it is higher for a larger mesh number. 

• The number of active grains per unit volume of media flow increases with percentage 
concentration of abrasives and reduction ratio, hence higher material removal and 
higher change in surface roughness are achieved. 

• Extrusion pressure, grain mesh size and hardness of workpiece material affect the depth 
of indentation of a grain in the workpiece material. Change in surface roughness value 
and material removal increase with extrusion pressure, number of cycles, reduction ratio 
and percentage concentration of abrasives. 

• A comparison of results predicted by theoretical model and simulation with RSA results 
reveals that results obtained from simulation are more close to RSA results than those 
obtained from the theoretical model. The original profile of the workpiece is taken to be 
Guassian in simulation, which is more close to the actual profile, than the triangular 
profile assumed in the theoretical model. 

Thermal aspects of AFM process 

• Specific energy for AFM process (10-110 J/mm^) is overlapping the range of specific 
energy for grinding (5-60 J/mm^). This justifies that machining action of AFM compares 
to grinding. 
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• Specific energy remains constant with abrasive mesh size. Us value increases with 
increase in hardness of the workpiece material and decreases with increase in applied 
normal pressure. 

• A simple model has been developed for heat transfer analysis in AFM process 

• It is found that the major portion of heat, developed at media-workpiece interface, 
enters into the workpiece. Hence, rise in temperature of workpiece is higher than rise in 
temperature of media. It is also observed that rise in temperature of workpiece 
increases with increase in extrusion pressure and number ot cycles. 

Application of neural network and genetic algorithm to AFM 

• AFM process is too complicated to warrant appropriate analytical model which are 
generally based on many assumptions. Hence, the possibility of using this neural 
network model for machined surface quality and material removal rate prediction for 
AFM process has been confirmed. 

• Appropriately trained network successfully synthesizes optimal input conditions for 
AFM process. 

• Genetic algorithm has been successfully applied for optimization of AFM process. The 
optimization results of the neural network coincide well with tlie results obtained by 
genetic algorithm. 
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7.2 Scope for Future Work 

In the present study, an attempt has been made to understand the mechanism of 

AFM process and its performance to a certain extent. Following are the areas in which 

further work can be carried out' 

1 As the present investigations are confined to cylindrical surfaces only, there is a need to 
study other shapes like flat and complex surfaces. 

2 AFM process is flexible in handling different types and shapes of tooling. Various types 
of tooling can be designed and experiments can be conducted for different shapes of 
workpieces especially for finishing complicated external shapes. 

3. There is a need to develop a mechanical device for efficient and uniform mixing of 
abrasive grains in the media. 

4. An extension of the theoretical model of media flow to three dimensional shapes and 
consideration of change in machining conditions of the AFM process with time are 
important issues that need to be investigated. Although the method has been defined 
for, and applied to, axisymmetric problems, the technique is easily extendible to three 
dimensional cases. 

5. Investigations in the area of media development are also .required The other material as 
carrier can be tested and developed economically. 

6. For better understanding of the process, investigations in the area of media viscosity 
and its flow characteristics are needed. 

7. As initial finish of the part is significant in AFM process, detailed study on this aspect 
will be of considerable interest for better understanding of the process. 

8. There is further scope to extend this work to finish workpieces of hard materials like 
ceramics and glass. 

9. The discussed neural network system is fairly general and can be extended to other 
abrasive processes to improve machining efficiency. The proposed off-line model can be 
applied for the on-line control of abrasive flow machining process. 
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10. Systematic modeling and analysis, optimum selection of process parameters, 
topography and micro-structural analysis of generated surface, abrasive grain wear 
phenomenon, and development of process control strategy are the main issues which 
need to be further investigated for continued developments and expended use of AFM 
process. 

11. The major constraint in the present work is the constant abrasive mesh size 
Experiments should be conducted in future, with different mesh sizes to understand its 
effect on the material removal and surface finish. The effect of back pressure on 
material removal and surface finish can also be studied 
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Experimental Measurements 


Table A.l Central composite rotatable design for 4 variables [25] 
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Table A-2 Material removal measurements 

, and ( indicate extrusion pressure in bar, number of cycles, reduction ratio and 
percentage concentration of abrasives in media respectively) 
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Table A.3 Surface roughness measurements 
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10 79 

8.02 

8.25 






2,17 


18,25 


13.17 


1.53 


1162 


8.34 







2.05 


16.63 


13.12 


139 


10.06 


7.65 







2.26 


17.92 


14 30 


1.47 


10,44 


9.02 


8 

35 

20 

93.75 

61 

1.65 

1.70 

14,73 

14.06 

9 . 46 . 

9.92 

0 . 90 . 


8,57 

8.72 

6 63 

6.95 






1.56 


12.30 


8 33 


0.9 


8.32 


5.00 







1.84 


15.41 


11.78 


1.17 


11.92 


9.10 







1.75 


13.83 


10.14 


0 87 


6.10 


7.10 



65 

01 




1.93 

14.00 

13.98 

10,10 

9.94 

091 

0 99 

6,89 

8.59 

4.85 

5.32 



n 





13.58 


9.74 


1.02 


8 72 


5.74 









13.51 


9 49 


1.04 


9.86 


5.11 




■ 





14.83 


10.43 


0.98 


8-79 


5.79 


10 

20 

15 

90 50 

66 

2.34 

2.31 

14.51 

13.93 

12.25 

12.13 

1.98 



12.62 

10 12 

10.59 






2.28 


14.56 


13,45 


2.08 




11.16 







2.37 


13.26 


11.25 


1.76 




9.66 







2.25 


13 42 


11.59 


2.13 




11.44 1 







































































































162 


Appendix A 






























































































164 


Appendix A 

















































































Appendix A 


165 


Table A.5 Active grain density for various concentrations of abrasives and 
mesh sizes 
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Table A.6 Measurements of temperatures at workpiece surface 


Ex 

No. 

Probe 

No 

Temperature at workpiece surface , "C 



Total 

temp 

rise, ”C 



O 

It 

A^=5 

N=10 

A'=15 

o 

il 

A/ =25 


1. 

1 

25 4 

25 8-25.9 

26.6 

27.1 



1 7o 


2. 

25 4 

25.8 

26.5 

27.1 



1 70 

2. 

1. 

25 9 

26.1-26 2 

26 4 




0.45 


2. 

25.9 

26.1-26 2 

26 3-26.4 




0.40 

3 

1 

25 3-25 4 

26.1-26.2 

26.9-27 0 




1 60 


2 

25 3-25.4 

26 2 

27 1 




1.80 

4 

1 

25 7-25 8 

26 1-26 2 

26.4-26.5 

26.8-26.9 

27.1-27.2 


1,30 


2 

25 7 

26 1-26.2 

26.3-26.4 

26.5-26.6 

26.8-26 9 


1.10 

5 

1 

26.0 

26.5-26 6 

26 . 9 - 21.3 

27.8-27.9 

28.4-28.5 


2.25 


2 

26.0 

26 5-26.7 

27 0-27.1 

27.8 

28.3 


2.10 

6 

1 

26.0 

26.3 

26.5-26.6 




0 55 


2 

26 0 

26 3 

26.5-26.6 




0.55 

IHH 


25.9-26.0 

26.7 

27.0 




1 05 

H 


25.9-26.0 

26.7 

27.0 




1 05 

8 

1 

26 2 

26.6-26.7 

26.8 

27.0 

27.2 


1,00 


2 

26.2 

26.6 

26.8 

27.0 



27,1 


0.90 

9 

1 

25.3-25.4 

26.0-26.1 

26.6 


27,5 


2.25 


2 

25 3-25.4 

26.1 

26.6 


27.5-27.6 


2.30 

10 


25.5 

"25?? 

25.8 




0.80 



25.5-25.6 

25.7-25.8 

1 

25.9 




0.80 

11 





29.4 



3.60 






29.4 



3.60 

im 


25.5 

26.1 

■■■[1 




0 60 



25.5 

26.1 





0.60 

13 

1 

25.6 

26.3 

26.9 

27.4 

27.8 

28.2 

2.60 


2 

25.6 

26.3 

26.8-26.9 

27.4 

27.8 

28.2 

2.60 

14 

1 


26.4 - 26.5 

26.9-27.0 

27.4-27.5 



1.45 


2 


26 4 

26.9 

27,3-27.4 



1.35 

15 

1 

25.8 

26.7-27.8 

26.8-26.9 

27.4 



1.60 


2 

25.8 

26.6 

26.8-26.9 

27.2 



1.40 

16 

1 

26.1 



g_|||||_ 



i.i) 


2 

26.2 






1.80 
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’Ex. Probe Temperature at workpiece surface , ^ C 
No. No. 


17 

1 

26 . 2 - 26.3 

26.9 


2 

26 3 - 26.4 

27 0 - 27.1 

18 

1 

26 . 1 - 26.2 

26 . 7 - 26.8 


2 

26 . 1 - 26.2 

26 . 9 - 27.0 

19 

1 

26 . 4 - 26.5 

27 . 2 - 27.3 


2 

26 . 4 - 26.5 

27 . 3-27 4 

20 

1 

26.1 

26.9 


2 

26.1 

26.9 

21 

1 

26.2 

26.9 


2 

26.2 

27.0 

22 

1 

26 1 

26.8 


2 

26.1 

26.8 

23 

1 

26 5 - 26.6 

26 . 7 - 26.8 


2 

26 . 6-26 7 

26 . 6-26 7 


1 

26.6 

27.4 


2 

26 5 

27.2 

25 

1 

26.1 

26 5 


2 

26.1 

26.5 

26 

1 

26.2 

27.1 


2 

26.2 

27.1 

27 

1 

26.1 

26 . 3-26 4 


2 

26.2 

26 . 4 - 26.5 

28 

1 

26.1 

27.3 


2 

26.1 

27.2 

29 

1 

26.1 

26.6 


2 

26.1 

26.6 

30 

1 

26.0 

27.1 


2 

26 0 

27.1 

31 

1 

26.1 

26.9 


2 

26.1 

26.9 


N=10 


27 . 4 - 27.5 

27 . 5 - 27 6 


27 . 3 - 27 4 

27 . 4 - 27.5 


27 . 8 - 27.9 

27 . 8 - 27.9 


27.5 

27.5 


27.4 

27.5 


27.4 

27.4 


27 . 4 - 27.5 

27 . 4 - 27.5 


28.2 

28.0 


rise, °C 


A^=15 


28 2 - 28.3 

28 . 3 - 28.4 


27 . 9-28 0 
28 . 0 - 28.1 


28 . 3 - 28.4 

28 . 4 - 28.5 


28.0 

28.0 


28.2 

28.3 


28.0 

28.0 


N=20 


iV =25 


27.6 

lie 


28 . 8 - 28.9 

28.8 
































































































Appendix B 

Viscosity of Media 


B.l Introduction 

Viscosity of media plays an important role in the performance of AFM process. At 
a particular velocity of piston, the viscosity of media affects the material removal and final 
surface roughness. The viscosity of media depends on type, mesh size and percentage 
concentration of abrasives, and type of carrier. From the literature survey, it is evident that 
there is a lack of investigation with regards to evaluation of rheological properties of the 
media. There is need to generate the data for viscosity of media of different types, mesh 
sizes and percentage concentrations of abrasives. The purpose of the present 
experimentation is to determine data for viscosity of media used in the AFM experiments 
reported in chapter 2 These viscosity experiments are necessary to calculate the material 
removal and surface roughness by theoretical model and simulation presented in chapters 
3 and 4. 

B.2 Capillary Viscometry 

In order to determine the rheological properties of the non-Newtonian fluids, many 
experimental methods have been developed. The viscometer developed based on the 
principle of capillary viscometry is the one used to determine the viscosity of the fluids in 
the present case. Hagan-Poiseuille studied capillary flow problems and derived relations 
between the flow rate and pressure drop considering the steady axial flow of an 
incompressible fluid in a circular tube [72]. 

In capillary viscometry (Fig.B.l), the media is forced through a fine-bore tube and 
the viscosity of the media is determined from the volumetric flow rate, applied pressure 
and tube dimensions. 

The following assumptions have been made in capillary viscometry 
• The flow is steady, incompressible, laminar, unidirectional and fully developed. 



Appendix B 


169 


• There is no slippage at the wall (this condition postulates that axial velocity, v = 0 at r 
= Aj , where r, is inner radius of the tube). 

• There are no external forces. 

• The tube is sufficiently long such that end effects are negligible. 


Load 


Rod-disc 

assembly 

Cylindrical tube 


Insulation 


Abrasive media 



Frame 


Fig. B.l Schematic diagram of a capillary viscometer 


For the fully developed laminar steady flow the expression for the two rheological 
functions, the shear stress and shear rate, are given by [72], 


21 , 

r = -— = f(r ) 

dr 


(Bl) 


(B.2) 


Where, r' is shear stress in N/m\ /, is length of the tube in m. Apis pressure drop 
across the capillary in N/m^ and u is axial velocity of the fluid in the tube in m/s 
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It is convenient to define the shear stress and shear rate at the wall of the tube 
(r = r,) in terms of measurable quantities 

. (B3) 


(B4) 


The flow of abrasive media starts after an initial shear stress is applied on it Hence, the 


expression for viscosity of media is given by 


(B.5) 


The initial shear stress of the material to initiate flow in the tube is found out 

by plotting the graph between the wall shear rate and wall shear stress (Fig.B.2) Using the 
method of least squares for fitting the data to a straight line, the initial stress is 
determined by extrapolating the straight line onto the x-axis. The intercept made on the x- 
axis at zero shear rate gives the value of . The slope of the straight line gives the 

reciprocal of viscosity. 



Fig. B.2 Relationship between wall shear stress and shear rate 
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B.3 Experimentation 

Media is filled into the tube of viscometer set-up upto its three-fourth length. 
Before filling the media into the tube, it is closed at the bottom by cylindrical cup. The 
media is then compacted by placing the 20 kg weight on the rod-disc assembly (Fig.Bl) 
which acts as a piston to force the media inside the cylindrical tube, for about half an hour. 
This procedure is repeated until the tube is filled with the media upto three-fourth of its 
length After properly filling the media in tube, the cylindrical cup is removed. 
Afterwards, the pressure is applied on the media by placing weights on the rod-disc 
assembly. Then, length of the extruded material is measured for a time interval of 2.5 
minutes. This procedure is repeated at different pressures. Table B. 1 shows the conditions 
under which the experiments have been conducted. 

Table B.l Experimental parameters used in calculation 
Radius of the tube = 47 mm 
time interval =2.5 min 
Temperature = 30 5 °C 


B.4 Calculation of viscosity 

Total Pressure gradient ( A p ) across the length of media is given by 
Ap=-^ (B.6) 

ra-t 

where, w is weight applied on the media which includes weight of rod-disc assembly. 

As extrusion starts, the pressure gradient changes across the length of the tube, 
because the length of the media filled in the tube initially, does not remain constant as 
soon as the load is applied. 

Assuming linear variation of /, with time as /,(t) = /„ - bt . (B.7) 

where, /,(t) = length of media in the tube after the time interval ‘t’, and b is a constant. 


As /, varies with time, is integrated over the interval 0 to t min. 
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1 r 

So, the actual value of 4 “ ~J (4 ” bX)dX 

^ 0 


(B.8) 


This actual value of 7,’ is used in the calculation of viscosity of media The 


constant b is determined by putting the measured value of /, at 0 and 2 5 min in eqn (B.7) 
If 4„,is extruded length of media in time^, then volumetric flow rate is given by 
following expression 

Q = (B.9) 

The values of shear stress and shear rate at wall at different pressures is determined 
by eqns (B 3) and (B 4) after substituting the values of A p and (3 • The initial stress 

is found out by plotting a graph between wall shear stress and shear rate as explained in 
sections B2 The viscosity of media is determined by eqn (B 5) after substituting values of 
and x\ 



Fig, B.3 Effect of percentage concentration of abrasives on viscosity of media 
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B.5 Results and Discussion 

The experiments are conducted on a capillary viscometer to determine the viscosity 
of the abrasive media at different concentrations It is evident from the Fig.B 3 that 
viscosity of media increases with increase in percentage concentration of abrasives. This 
can be explained as follows' As the % concentration of the abrasives in the media 
increases, it decreases mobility of the particles migrating from the tube wall Due to this, 
the amount of material extruded from the tube decreases with increase in percentage 
concentration of the abrasives in the media, which ultimately decreases volumetric flow 
rate and the wall shear. There is a decrease in shear rate due to decrease in the volumetric 
flow rate Due to this reason, there is an increase in viscosity with increase in percentage 
concentration of abrasives in the media. 


B6. Experimental Data 


Abrasive concentration = 56 % 
Room temperature = 30 5 °C 
Initial stress = 13 74 KN/m^ 


s 

No 

Load 

Applied 

(kg) 

Extrude 

d length 

m. 

Average 

length of 

the media 

r(m) 

Wall shear 

stress 

(KN/m^) 

i 

Shear 

rate 

^ (s‘) 

Viscosity 

A 

(KN-s/m^) 

Av. 

Viscosity 

(KN-s/m^) 


53 7 

0 0250 

0.4705 

39.37 

0.1418 

182.72 

185 40 

■ 

58 7 

0.0275 

0 4670 

43.49 

0.1560 

190.68 


3. 



HH 

47.30 

0 1787 

187.77 


4 

68,7 

0.0360 

0,4650 


BBi 

185.41 
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Abrasive concentration = 61 % 
Room temperature = 30 5 °C 
Initial stress = 14 78 KN/m^ 


s 

No, 

Load 

Applied 

(kg). 

Extruded 

length 

m. 

Average 

length of 

the media 

r(m) 

Wall shear 

stress t\ 

(KN/m^) 

Shear 

rate 

r (s'*) 

Viscosity 

(KN-s/m^) 

Av 

Viscosity 

(KN-s/m^) 

1- 

53.7 

0 0205 

0 4735 

36.92 

0 1163 

207.61 

205.60 

2 

58 7 

0.0240 

0 4705 

42 99 

0.1361 

207.19 


3 

63 7 

0 0280 

0 4680 

47.11 

0.1588 

203,55 


'4 

68 7 

0,0315 

0 4640 

51.46 

0.1787 

205 28 



Abrasive concentration = 66 % 
Room temperature = 30.0 °C 
Initial stress = 15.50 KNW 


S. 

No. 

Load 

Applied 

(kg). 

Extrude 

d length 

m. 

Average 

length of 

the media 

I(m) 

Wall shear 

stress r\ 

(KN/m^) 

Shear 

rate 

r (s'*) 

Viscosity 

P- 

(KN-s/m^) 

Av. 

Viscosity 

(KN-s/m^) 

1. 

53.7 

0.0175 

0 4760 

38.59 

0 0992 

232.59 

231.87 

1 - 

2. 

58,7 

0.0200 

0.4740 

42.48 

0.1134 

237.79 

3, 

63.7 

0.0240 

0,4720 

46.50 

0.1362 

227.68 

4. 

68 7 

0 0270 

0.4690 

50.67 

0.1532 

229.44 
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Abrasive concentration = 71 % 
Room temperature = 30.5 °C 
Initial shear stress = 16.78 KN/m^. 


s. 

No. 

Load 

Applied 

(kg). 

Extrude 

d length 

m. 

Average 

length of 

the media 

r(m) 

Wall shear 

stress 

(ICN/m^) 

Shear 

rate 

r (s') 

Viscosity 

(KN-s/m^) 

Av 

Viscosity 

(KN-sW) 

1. 

53.7 

0.0140 

0.4685 

39.07 

0.0794 

280.68 

273.35 

2. 

58.7 

0.0165 

0.4750 

42.23 

0.0936 

i 

271.89 


3. 

63.7 

0.0195 

0.4720 

46.28 

0.1106 

266.61 


4. 

68.7 

0.0215 

0.4700 

50.23 

0.1219 




Abrasive concentration = 76 % 
Room temperature = 30.0 °C 
Initial shear stress = 17.66 KN/m^ 


S 

No. 

Load 

Applied 

(kg). 

Extrude 

d length 

m. 

Average 

length of 

the media 

E(m) 

Wall shear 

stress 

(KN/m^) 

Shear 

rate 

r (s *) 

Viscosity 

M 

(KN-sW) 

Av. 

Viscosity 

(KN-s/m^) 

1 

53.7 

0.0115 

0.4740 

38.51 

0.0652 

319.58 

308.07 

2. 

58.7 

0.0145 

0.4710 

42.50 

0.0822 

302.01 


3. 



0.4700 

46.32 

0.0936 

306 22 


■1 

68.7 

0.0190 

0.4665 

50.48 

0.1078 

304.47 
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Appendix C 
Specifications 


Workpiece 

Material = mild steel {025C),H^ = 2177.80 N/mm^, 

38.5 mm, 43.40 J/kg-K, k^=6Q.5 W/m-K, 

7854.0 kg/m^ 

Abrasive [40] 

Type = SiC, M = 50-60, 55, 1400.0 J/kg-K, 

= 100.0 W/m-K, = 3220 0 kg/m^ 

Carrier - Polyborosiloxane 

[22] 

= 434.0 J/kg-K, k= 0.22 W/m-K, p^= 985.0 kg/m^ 

Media Cylinder 

/'^= 43.75 mm, 1= 80.0 mm 

Assumed 

/=0.2[41] 

value of coefficient of friction 

A = 0.002 s [33] 

Relaxation time 
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Appendix D 

Data from AFM experiments [19] 


SN. 


C % 

M 

N 

MRR 

Ra 

1 

40.60 

45.0 

100 0 

50 

1.464 

1 63 

2 

51 50 

45 0 

100.0 

50 

2.270 

1.08 

3 

65 20 

45 0 

100 0 

50 

3 260 

0.70 

4 

81.20 

45 0 

100 0 

50 

5.684 

0.47 

5 

51 50 

45 0 

150 0 

50 

1.957 

1 14 

6 

51 50 

45 0 

180 0 

50 

1.442 

1.18 

7 

51.50 

45.0 

240.0 

50 

1.236 

1.23 

8 

40.60 

41.0 

150.0 

80 

0.869 

1.80 

9 

65.20 

41 0 

150.0 

40 

2.147 

0.88 

10 

51.50 

33 0 

100.0 

20 

1 545 

1.94 

‘■ni 

51.50 

33.0 

100.0 

40 

0 901 

1.54 

12 

51.50 

33 0 

100 0 

60 

0.858 

1.52 

13 

51.50 

33.0 

100.0 

80 

0 837 

1.48 

14 

51 50 

33,0 

100.0 

100 

0.772 

1.42 

15 

40.60 

45 0 

150.0 

80 

1.149 

1.59 

16 

65 20 

35.0 

180.0 

40 

1 111 

1.12 

17 

81 20 

35 0 

240 0 

60 

1.119 

0.71 

18 

51.50 

37 0 

100.0 

20 

1.802 

1.75 

19 

51.50 

37 0 

100.0 

40 

1.545 

1.44 

20 

51.50 

37.0 

100.0 

60 

1.459 

1.48 

21 

51.50 

37.0 

100.0 

80 

1.481 

1.33 


51.50 

37 0 

100.0 

100 

1.442 

1.25 

23 

40.60 

45 0 

150.0 

100 

1.102 

1.51 

24 

65.20 

33 0 

180.0 

80 

0 814 

1.03 

25 

81.20 

33 0 

240.0 

40 

1.017 

0 84 

26 

51 50 

41 0 

100.0 

40 

2.704 

1 34 

27 

51.50 

41.0 

100.0 

60 

2 661 

1.15 

28 

51.50 

41.0 

100.0 

80 

2511 

1.09 

29 

51 50 

41.0 

100.0 

100 

2.214 

1 04 

30 

65.20 

35 0 

150.0 

60 

1.222 

1.01 

31 

81 20 

41.0 

180.0 

20 

2 692 

0.72 

32 

51.50 

45.0 

100.0 

40 

2 832 

1.28 

33 

51.50 

45 0 

100.0 

60 

2 575 

1.03 

34 

51.50 

45 0 

100.0 

80 

2.446 

0.95 

35 

51.50 

45 0 

100 0 

100 

2.369 

0.87 


where, Vp is piston velocity in cni/min, M is mesh size of abrasive grains, C is 
percentage concentration of abrasives by weight, N is number of cycles, MRR is material 
removal rate in g/minx lO"^, and is surface roughness in /an . 
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